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SUMMARY 


A comprehensive analytical design procedure for the installation of simple 
pendulums on the blades of a helicopter rotor to suppress the root reactions is 
presented. To achieve this goal, a frequency response analysis is conducted of typical 
rotor blades excited by a harmonic variation of spanwise airload distributions as well as 
a concentrated load at the tip. 

The structural modeling of the blade includes elastic degrees of freedom in flap and 
lead-lag bending plus torsion with a hingeless hub constraint. Simple pendulum 
absorbers are individually considered for both flap and lead-lag types of motion. The 
inertial reaction forces and moments at the pendulum hinge are formulated. On the 
basis of a rational ordering scheme the general nonlinear equations of motion for the 
rotor-pendulum system are linearized in the perturbation elastic blade displacements 
and pendulum angle. These linearized equations include the effects of spanwise airload 
distributions associated with the elastic deformations and the cyclic pitch environment 
of forward flight. A quasi-steady aerodynamic representation is utilized in the 
formation of these airloads. The solution of the system equations is based on their 
representation as a transfer matrix. 

The numerical results presented here pertain to both uniform and nonuniform 
hingeless rotor blades. These results include the effect of pendulum tuning on the 
minimization of the hub reactions. It is found that a properly designed flapping 
pendulum attenuates the root out-of-plane force and moment whereas the optimum 
designed lead-lag pendulum attenuates the root in-plane reactions. A properly tuned 
pendulum can attenuate the vibratory loads by generating appropriate forces at its 
attachment point with the blade. These for :es redistribute the loads on the blade so 
that only a small portion of the reactions are transmitted to the hub. 

For optimum pendulum tuning the parameters to be determined are the pendulum 
uncoupled natural frequency, the pendulum spanwise location and its mass. It is found 
that the optimum pendulum frequency is in the vicinity of the excitation frequency. In 
the case of an off-tuned pendulum the result can be either a slight-to-moderate 
degradation in pendulum absorber effectiveness or a drastic increase in hub reactions. 
Although the uncoupled natural frequency of the pendulum is independent of the mass of 
the pendulum, a proper choice of the mass is required to generate an optimum force to 
attenuate the root reactions. 

A pendulum can be tuned and its optimum mass determined by excitation with a 
concentrated simple narmonic load at the tip. However, it is necessary to utilize 
distributed airloads to accurately determine the attenuation of the root reactions. The 
damping at the pendulum hinge has a small effect on the hub reactions once the 
optimum pendulum tuning is established. For the optimum pendulum a parametric study 
is conducted. The parameters varied include prepitch, pretwist, precone and pendulum 
hinge offset. 


INTRODUCTION 

Vibration has always been a significant detriment to helicopters which is not 
surprising in view of the fact that a major percentage of a helicopter's mass consists of 
rotating components such as engines, rotors, gearboxes and drives. These revolve at 
different frequencies, though most are multiples or submuitiples of the main rotor 
frequency, and any out-of-balance forces are easily transmitted to the light airframe. 
There is also the possibility of damaging beat frequencies arising in those few cases 
where rotational rates are not synchronous. 


In addition to the inertially induced vibrations, the unsteady air loads on the main 
and auxiliary rotors significantly contribute to the vibratory environment, These 
unsteady aerodynamic reactions are caused by the local periodic dynamic pressure on 
the main rotor blades in forward flight and the cyclic pitch changes imposed by the 
rotor head. In response to these unsteady loads the elastic blades are continually 
deforming which induces additional unsteady aerodynamic loading. 

An increasing demand for the reduction of helicopter vibration has been dictated by 
the expansion of flight envelopes coupled with more stringent requirements for crew 
and passenger comfort as well as improved reliability and maintainability. 

The rotor system, which transmits the vibratory airloads to the fuselage through 
the rotor shaft, is one of the most significant contributors to the vibrations of the 
helicopter. The rotating blades create vibratory airloads containing ail harmonics of 
the rotor rotational frequency. As the loads from the individual blades combine at the 
rotor hub, some harmonics cancel each other while others are additive. The loads that 
are additive are passed from the blades to the pylon and then to the cabin through 
complicated dynamic load paths. These loads are felt as vibratory forces and moments 
whose frequencies are integer multiples of the blade passage frequency (number of 
blades times rotational frequency). 

Dynamic systems that reduce the vibration in helicopters may be classified into 
five groups: 

(1) Excitation reducers (blade aerodynamics) 

(2) Attenuators (blade dynamics with or 

without hub motion) 


(3) Absorbers (hub dynamics) 
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(5) Cabin Suppressors (fuselage dynamics) 

These systems are described below and specific examples are given. 

Excitation reducers are devices that either reduce aerodynamic forces on the 
blades, or alter the hub moments and shears by generating counter-influencing 
aerodynamic forces on the blades, Higher harmonic blade pitch control devices (ref. 1) 
fall into this category. In reference 1 a wind tunnel test has been used to evaluate 
higher harmonic blade pitch for the reduction of helicopter vibration. The higher 
harmonic pitch was obtained by higher harmonic oscillation of the swashplate. The test 
results indicated reduction of fuselage vibratory loads. 

The attenuators are those dynamic devices that result in low hub shears and 
moments while the blades are subjected to specified aerodynamic excitation forces. An 
example of this category is a midspan or tip weight that may move the natural 
frequency of a mode away from exciting harmonics or reduce the modal forces at the 
blade root. 

Absorbers are devices that are incorporated at the hub (or blade root) and absorb a 
large portion of the hub shears or moments causing transmittal of only the remaining 
forces to the mast. Examples are the simple and bifilar pendulums. 

Isolators are generally integral to the pylon assembly and reduce greatly the forces 
and moments transmitted to the fuselage with respect to the forces applied on top of 
* the mast. 


Cabin suppressors reduce the vibration in the cabin for prescribed forces and 
moments applied by the pylon on the fuselage. 

One of the five categories, the pendulum absorber, has proven to be quite 
successful as reported in references 2 to 6. Such a technique consists of mounting a 
simple or bifilar pendulum on each rotor blade near the hub. The blade-mounted 
pendulum absorbers modify the response dynamics of the blades, and reduce the 
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helicopter vibrational levels. All reported work has consisted of an expensive and time- 
consuming flight test program to establish a set of acceptable pendulum parameters. 
As the flight tests become more and more expensive, an analytical study is needed for 
pendulum design. 

Murthy and Hammond (ref. 7) have studied the effect of a pendulum absorber on 
the natural frequencies and mode shapes of a rotor blade. A vibration analysis of rotor 
blades undergoing coupled flapwise bending, chordwise bending, and torsional vibrations 
with a spherical pendulum absorber was presented. A portion of the investigation 
included a reduction of the spherical pendulum to a simple flapping pendulum. It was 
shown that the pendulum absorber displaces original bladt natural frequencies that are 
in the vicinity of the pendulum tuning frequency(ies). The spherical pendulum 
introduces two modes between the displaced original modes of the blade, whereas the 
simple flapping pendulum introduces one mode. Mode shape results were not presented 
in the report. It is worth mentioning that the equations of motion were based on the 
linear analysis presented by Hpuboit and Brooks (ref. 8), and no damping terms were 
included. Also the damping at the pendulum hinge was not considered. 

The principal aim of the present study is to establish comprehensive analytical 
design procedures and criteria for the installation of simple pendulums on the blades of 
a helicopter rotor to suppress the force and moment reactions at the hub during forward 
flight conditions. 

The objectives of tha present investigation can be stated as follows: 

(i) Develop a mathematical model to represent the blade-pendulum system. A 
single nonuniform rotor blade with a hingeless hub restraint undergoing coupled flapwise 


bending, eherdwise bending, and torsional vibrations is considered. Simple pendulum 


absorbers are individually treated for both flap and lead-lag types of motion. The blade 
is excited by an azimuthal harmonic variation of spanwise airload distributions 
associated with the elastic deformations and the cyclic pitch environment of forward 
flight, 

(2) Find the dynamic response characteristics of the blade-pendulum system, using 
the transfer matrix method (ref, 9). 

(3) Determine the optimum pendulum tuning to suppress the hub reactions. This 
entails the minimization of these reactions by appropriate variations of the pendulum 
parameters for a given excitation frequency. The pendulum parameters include the 
uncoupled pendulum frequency, hinge spanwise loaction and pendulum mass. 

( 4 ) Conduct a parametric study of the optimum tuned configuration. The 
parameters to be varied include pendulum hinge offset, precone, prepitch and pretwist. 
The intention of this investigation is to document the effects of these parameters on 
the optimum configuration previously established and thereby provide useful design 
criteria for future installations of pendulum absorbers. 


SYMBOLS 
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Ct 1 


airfoil lift-curve slope 

reference frame which rotates with speed ft with 
respect to the stationary inertial frame Off’ 
pendulum hinge friction coefficient, equation (35) 
rotor thrust coefficient 
blade chord 
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airfoil profile drag coefficient 

Young's modulus of elasticity 

mass centroid offset from elastic axis, positive when 

forward of the elastic axis 

aerodynamic loads vector, equation (133) 

components of aerodynamic forces per unit length in y', <\ 

y', z' directions 
shear modulus 
acceleration due to gravity 

"fixed" unit vectors with & in direction of rotation of * 

01n &? , figure A-l 

blade cross-section moment of inertia about y' axis 
blade cross-section moment of inertia about z' axis 
unit vectors assolcated with undeformed blade coordinate 
system, figure A-l 

unit vectors associated with deformed blade coordinate 

system, figures 1, A-2 

torsional rigidity constant 

radius of gyration of blade cross section 

mass radius of gyration of blade cross section 

mass radii of gyration about ri and £ axis, respectively 

components of generalized nonconservative (aerodynamic) 

forces in y and z directions 

length of pendulum 

mass of the pendulum 

twisting moment about the x and x' axes 

bending moment about y and y' axes 

bending moment about z and z' axes 

generalized nonconservative (aerodynamic) moment about 

x axis 

mass per unit length of the blade 

time independent average normal force on the pendulum 

bearing, equation (35) 

fixed point in^and^at the root of the blade, figure A-l 
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W 

X, Y, Z 


x, y, z 


x', y', z' 


{Z(x,t)} 

a s 

^pc 

e 
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point transfer matrix of the pendulum 

length of blade 

inertial reference frame 

transfer matrix of the blade from 0 to x 

initial radial blade tension, equation (10) 

time 

radial, tangential, and perpendicular components of 
velocity for blade airfoil section 
elastic displacements in the x,y,z directions, 

respectively, figure A-2 
wind velocity 

axial force (tension) in the x, x’ directions, 
respectively 

shear force in the y, y' directions, respectively 
shear force in the z, z' directions, respectively 
induced downwash velocity at rotor, positive 
downward 
pendulum weight 

inertial axis system with origin at hub centerline and 
Z normal to hub plane, figure A-l 
mutually perpendicular axis system with x along the 
undeformed blade and y toward the leading edge, 
figures A-l 

a set of Cartesian coordinates with origin at a point 
along the deformed blade with x' remaining tangent 
to the elastic axis and y’ and z' along principal axes 
for the cross section, figures 1, A-2 
state vector, equation (1) 
shaft angle 

precone angle, figure A-l 

small parameter of the order of the bending slopes, 
also airfoil section pitch angle with respect to free- 
stream velocity, figure 7 

sectional coordinate normal to n axis at elastic axis 
(same as z'), figure 1, also damping ratio 
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Subscripts: 

A 

P 

Superscripts: 

L 

R 

T 

-1 

Notations: 

( )' 


( * ) 
<“) 


sectional coordinate corresponding to major principal 
axis for a given point on the elastic axis (same as y'), 
figure 1 

pitch angle, equation (125) 
collective pitch, equation(125) 
cyclic pitch components 
pretwist, equation (126) 

pretwist per unit span, positive when tip angle is 
larger, equation (126) 
inflow ratio, positive upward 

slope of deflection curve in the plane of rotation, 
dv/dx 

advance ratio, jj^SJR = Vj cos a s 

slope of deflection curve normal to the plane of 

rotation, dw/dx 

mass density of blade, also mass density of air 

elastic twist about the elastic axis 

blade azimuth angle measured from downwind position 

in direction of rotation 

rotational speed of the blade 

uncoupled natural frequency of the simple pendulum 


at the pendulum hinge location 
at the pendulum mass location 

immediate left 
immediate right 
transpose of matrix 
inverse of matrix 

differentiation with respect to x, also deformed 
coordinate system 
differentiation with respect to time 
amplitude of simple harmonic quantity 
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BLADE EQUATIONS OF MOTION 


The equations of motion which are used to represent the structural dynamic 
characteristics of the helicopter blade are based on the development reported by 
Hodges and Dowell In reference 10. Tne mathematical model chosen to represent the 
rotor blade consists of a straight, slender, variably twisted, nonuniform beam i h'eh can 
undergo combined flapwise bending, lead-lag bending and torsion. The elastic axis? the 
mass axis, and the tension axis (area-centroid axis) are taken to be noncoincident. The 
elastic axis is inclined to the plane of rotation at a small angle in order to accommodate 
any built-in coning (precone). The use of an actual helicopter rotor blade as a model 
would be a more formidable task because of the nonhomogeneous, anisotropic 
structures found in typical blades. 

Several orthogonal coordinate systems and coordinate transformations (Appendix A) 
are employed In the derivation of the equations of motion. The beam cross section is 
shown in figure 1 before and after deformation. The n and 5 axes are the principal axes 
of the cross section for the shear center. The cross section is assumed to be symmetric 
with respect to the r\ axis. The deformed beam is shown In figure 2(a) with force and 
moment resultants acting on the face of a cross section. The x, y, z coordinate system 
is the preconed undeformed system. At any point along the deformed beam, x' is 
tangent to the deformed elastic axis. The y' -and z' are identical to the n and £ axes, 
respectively, when the beam is deformed. 

A mathematical ordering scheme which is consistent with the assumption of a 
slender beam is adopted for the purpose of systematically discarding elastic and 
dynamic terms which are of higher order in the equations of motion. The ordering 
scheme is consistent with the small deformation approximation in which the strains are 
negligible compared to unity and the dimensionless axial deflection u/R is generally 
taken to be of the same order of magnitude as the square of v/R or w/R and thus is 
small with respect to unity. These assumptions can be systemized by introducing e , a 
parameter of order v/R or w/R. The ordering scheme associated with the above 
consideration is summarized in Table 1. Within the force and moment expressions, 
terms of order e 1 are neglected with respect to umty. Thus, if the largest terms of the 
force expression are 0( e 7 , then all terms of 0( e 7 are retained (first-order terrra), all 
terms of 0 ( e p ) are retained (second-order terms), and generally terms of 0 ( e 7 are 
discarded. 

The blade equations of motion will be written in terms of the state vector { Z } 
where 


{Z(x,t)} T = [vw * P V M x M y M z V x V y V Z J (1) 

It may be noted that the above state vector does not Include the axial deflection, u. In 
order to obtain a set of linear equations for the state vector components, u must be 
eliminated. The components of the state vector, { z} , can be chosen several ways, but 
they are chosen here such that they represent the physical quantities of deflections, 
slopes, moments, and shears. This is not absolutely required but highly preferable for 
application of the transfer matrix to obtain the dynamic response characteristics. 

The equilibrium of the forces and moments that act on a differential beam element 
is now considered. In this consideration the element is formed by slices parallel to the 
yz-plane, because this choice leads to rather simple results. The forces that act on such 
an element are shown in figure 2(a); the moments, in figure 2(b). The quantities p , p , 
p , q , q and q are resviltant force and moment loadings. Summation of the moments 
aBouf the x-, and z-directions lead to the following equilibrium conditions for 
moment and shear: 


7 


JOK JO!* -e- J® 1< 5° IC 


TABLE 1. - ORDERING SCHEME FOR ROTOR 
BLADE PARAMETERS 


0( e 2 ) 

0 s 

0(1) 

0( e) 

^ pc = 

0(e) 

* 

0( e ) 

u 

,51* 

o ( G ) 

0( c) 

k^ 

m 

R s 

0(g) 

0(1 ) 

i 

R( )' = 

0(1) 

0(g) 

S ( '> = 

0(1) 

0(g) 
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M' x + V 2 v* -Vy w* + q x 

a 0 

( 2a) 

M'y ~V 2 fV x w' + q y 

= 0 

( 2b) 

M 'z f V y -V x V' + q z 

a 0 

( 2c) 

V 'x + Px 

a 0 

( 2d) 

W 

= 0 

( 2e) 

V'z + Pz 

a 0 

( 2f) 


The applied loads p , p , p , q , q , q due to inertial and aerodynamic loadings a<e 
functions of v, w, (J) ancra longuudinarradial deformation u in the x direction. The only 
effective aerodynamic loading to be considered in this analysis will be forces in the y 
and 2 directions given by L and L , respectively, and moment about the x - axis given 
by M . . The inertia loadings, which consist of the forces and moments which oppose 
accelerations of the blade elements, are derived in reference 10 whereas the 
aerodynamic loadings L , L , and are left in this symbolic form. These 
aerodynamic loads are derived later in this text. The procedure that follows is 
employed in the derivation of the inertia loading which is given in reference 10. The 
acceleration of any mass particle on the vibrating, rotating blade Is derived and the 
components in the x-, y- and 2 - directions are obtained. These component accelerations 
include terms for the Coriolis, transverse and centripetal accelerations along with 
higher order terms. The inertia force and moment loadings are then obtained by 
integrating over the cross section. By dropping higher order terms in accordance with 
the adopted ordering scheme, the following resultant loadings to second order are 
obtained: 


% = 


- m( fi 2 x + 2 ft V)-m R w 

( 3a) 

a -mV + me <j> sin 0 + m fi 2 [v + e cos (0 + <}>)] 

( 3b) 

+ 2m fi( $ „w-Ci) + 2m fie (v' cos 0 + w'sin 0 ) + L 
pc V 


= -mw -me IjJcos 0 -m fi $ pc x -2 mfi$ pc v + L w 

(3 c) 

• . 2 2 
a me(v sin 0 -w cos 0 ) -m fi ev sin 0 -m fi e $ x cos 0 

pc 

( 3d) 

~mk^(j) m fi 2 (k^ -k^ ^ ) sin 2( 0 + <f> ) + 


2 

= me( fi x sin( 0 + <J>) + 2 fiv sin 0 ) 

( 3 e ) 

2 

= -me( fi x cos(0 +<()) + 2 fiv cos 0 ) 

(30 


The underlined higher-order terms could be dropped consistent with the ordering 
scheme, but this would eliminate torsional inertia. In the special case of rotor blade 
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configurations of very low torsional rigidity, these terms may contribute substantially 
to the magnitude of the torsiona' natural frequency? lienee, they will be retained (ref. 
10 ). 


w 


There are now 6 equilibrium equations, equations (2a) - (2f), and 10 unknowns, u, v, 
w, <J>, M , M , M , V , V , and V . In addition, 4 force and moment deformation 
relations are requited to express M f M , M , V In terms of u, v, w, <t> ♦ The force and 
moment resultants are expressed* in y tcrms of the deformations by resolving the 
distributed stresses into a resultant force and moment system acting at the elastic axis. 
These relations given in reference 10 are written in the deformed body axis system as 


i v' 2 w' 2 

V* - E M u ' + T + -2 


* k A 0 ' 4 >’-c a |y" cos (0 + <>) + w" sin (0 + $ )]| 


( *a) 


M x' 


2 2 

GO <j>» + EAk 2 (0 + $)' (u' + j- + ^~) + EB*^9 ,2 (j>» 


( 4b) 


EB| 0'(v" cos 0 + w" sin 0 ) 


My, = EI yt [v" sin ( 0 + <t>) -w" cos (0 + c|>)] (4c) 

2 2 

M z , « EI z , [v" cos (0 + <},) + w” sin (0 + * )] ~EAe A (*+$**%) -EBJ 0 f (4d) 


where B| and B| are section constants defined by 


B i 5 ff s 2)2dr i d s 

A 

B 2 n 2 +? 2 )dnd5 

A 

and A is the blade cross-section area effective in carrying tension. 

The underlined terms may be neglected according to the ordering scheme. 
However, as discussed before, these terms are important for configurations with low 
torsion stiffness. In equations (4a) - (4d), the terms which come from the warp effect 
are omitted due to its relative smallness for closed cross sections. 

The force and moment deformation relations, equations (4a) - (4d), can be readily 
transformed into the undeformed body axis system using the transformation matrix of 
Appendix B. The resultant moments are; 
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,2 


.2 /„ . . v' 2 . w'. 


M x = GJ (J)’ + EAk* (0 + <j>)‘ (u' + y + 


+ EBJ 0' 2 t})' -EBg 0* (v" cos 0 + w" sin 0) 


,2 ,2 

- EAo A (u* + y* + ^-) (v* sin 0 -w' cos 0 ) 


+ (EI Z , -Ely,) (v V' -w'w") sin 2 0 + (v'w" 


w'v") sin 2 0 J + Ely, v'w" -EI Z , v'V 


v' 2 . w' 2 


M = GJ (j>V + EB| 0 ' <j>' sin 0 + EAe^ (u 1 + -g- + ^*) 


• sin (0 + <j)) -(EI Z , -Ely,) [w- sin 2 <0 + <j)) 
+ ~v" sin 2(0 + (),)] -El , w" 


M. 


= GJcJj'w' -EB| 0' (J)' cos 0 -EAe A (u' 


2 2 

+ T + T>cos(0 + <|») + (EI 2 ,-EI ,) 


• j [w" sin 2 ( 0 + <}> ) -v" sin 2 (0 + 4>)J + EI z ,v" 


The resuitant axiai force in the undeformed body axis system is 


V x = V x , -V , jjv* cos ( 0+ <j>) + w' sin ( 0 +<j> )J -V 2 , [-v' sin (0 +«j> ) + w' cos (0 + 


( 5a) 


( 5b) 


( 5c) 


+>] 
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For rotor blade applications the second and third terms on the rifht-hand side are 
negligible to second order (ref. 10)*, thus*, 

V x st V x , = EA ju'+y + y- 2+,< a 0 ' ( i ,, “ e A [V’cos (0 + <J>) + w" sin (0 | (6) 


In the above expressions all terms above second order have been neglected. As a final 
step the following definitions are introduced. 

v* a „ (7a) 


w' 5 V (7b) 

Equations (2), (5), (6), and (7) represent 12 independent equations for the 11 state 
variables and the elastic elongation, u. In order to obtain a set of linear equations for 
the state vector components, u must be eliminated. This will be accomplished by 
defining the time independent tension, T , as 

s 

T 0 = / m fi 2 x d x 

In the process of linearization, the following approximations will be made 


sin (0 + $) ss sin 0 + cos 0 


cos ( 0 -h <b ) a cos 0 - sin 0 

2 sin 2(0 + tt>) a ^ sin 2 0 + 4 > cos 2 0 

2 2 
sin (0 + 4 ) s sin 0 + <j> sin 2 0 


( 8 ) 


Equation (6) can be rewritten as 
,.2 „.»2 V 


, V w 

u + 2 + 2 " EA _, 'A 


y 2 

- I< A 0 1 (j) ' + e A y' cos 0 + e^v 1 sin 0 


( 9) 


The left-hand side of equation (9) appears in equations (5) with both constant and linear 
coefficients. When the coefficient is constant it shall be represented as written above 
but with a linear coefficient it will be approximated as 


(u* + 



W |2 v 

2 'o 


C* 



< 10) 
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In the process of linearization, all steady-state values are neglected except the 
initial radial tension, and only the linear terms in the state variables are retained. 
These steady-state values would correspond to an equilibrium flight condition in hover. 
It is anticipated that these equilibrium hover values would have little effect on the 
dynamic response characteristics in the forward flight condition of this study. The set 
of linearized first-order differential equations are given below 


V * w (11) 

w' cs v (12) 

<j>' (GO + k 2 T o + EB* 0 ,2 ) s -!< 2 a 0 ' V x + EA(k 2 0 ') 2 <j> * 

-EAk^ 0 'e A ( y' cos 0 + v'sin 0) + EBJ 0 '( y'cos 0 
+ v'sin 0) + w ysinO - vcos 0) + M x (13) 

y ' (EL, cos 2 0 + E! , sin 2 0 ) = M + EB* 0 ' d>' cos 0 + e , cosO V 
z y z /» a x 

- EAk A e A 0 * <j)' cos 0 + EAe A ( y'cos 2 0 + v ' sin 0 cos 0) 
-e A T q 4> sin 0 -|<EI 2 , -EI yt ) V ’ sin 2 0 ( H ) 

V ' (EI y , cos 2 0 + EI 2 , sin 2 0 ) = -M y + EB* 0 ' <}> ' sin 0 + e A V x sin 0 

-EAI< a e A 0 ' (j)' sin 0 + EAe A ( y'sin 0 cos 0 + v' sin 2 0 ) 

+ e A T 0 cos 0 " k < EI z' “Ely,) y ' sin 2 0 (15) 


M' = 


' 2 . 


mfl e x( v cos 0 - y sin 0 ) + (e A T q cos 0 )'v -(e A T Q sin 0 )' y 
+ mk^ $ + me [( fi 2 v -v) sin 0 + w cos 0 ] 

( 16 ) 


ip 


M' 


V„ -T v -2 mil ev sin 0 -mfi 2 ex <b cos 0 
z o T 


( 17) 


13 


( 18) 


2 

Mlj, a T Q y-V y + 2 m ftev cos 0 -mft ex <j> sin 6 

\jj = -2 mftv (19) 

2 2 

V'y = mV -me $ sin 0 -m ft v + m ft e <j> sin 0 -2 m ftB pc w 

i 

-2 m fie jjcos 0 -2 mfte vsin8-L v (2 0) 

V s m\V + me $ cos 0 + 2 m fiB v -L (21) 

/* pQ W 


It may be noted that the three equations for the elastic slopes (eqs. (13), (14), and 
(15)) are coupled in their derivatives and must be solved simultaneously, These 
equations are decoupled by a straight forward application of Cramers rule to obtain a 
set of linear first order differential equations in terms of the state vector so that the 
transfer matrix method can be applied as discussed later in this study. 

Excitation of the biade is due to various simple harmonic airloads. The solution to 
the set of equations will be treated as individual contributions associated with each 
harmonic of the excitation loading. Such a solution for the nth harmonic of the state 
vector can be written in complex form as 


{ Z n (x,t) } » Re((Z n (x)J e lut ) (22) 

where { Z n (x)} is the amplitude vector of the nth simple harmonic response, and 

a) = n ft = the frequency of the excitation loading 
i =VT 

The above solution will be discussed in more detail with respect to the aerodynamic 
loading. Substituting equation (22) into equations (11) to (21) yields the following matrix 
equation for { Z} , the nth harmonic, 


{ Z (x) } ' = [A (x)] { Z (x) } + { F q (x)} 


where { F } is the amplitude vector of the aerodynamic loads, defined by 


{ F } T = I 0 0 0 0 0 ~M , 0 0 0 -L - L I 

1 a 1 l. cf> v w-J 


14 


( 23 ) 


Note that the matrix, [a], Includes complex terms due to retaining the linear rate 
terms (Coriolis intertial forces). 

Evaluation of the order of magnitude of the individual elements of the complex 
matrix [,A] is made. Consistent with the adopted ordering scheme, only terms of first 
and second order are retained. It may be noted that in the process of decoupling eqs. 
(13), (14), and (15), terms of higher order are encountered. The resulting nonzero 
elements of [A] become 


A 14 


1 

A 25 

= 

1 

A 34 

= 

El .(EI^-EAe* ) e A T o sin 0 /D 

A 35 

S 

-A 3 ^ cot 0 

A 36 

= 

EI y ,(EI z ,-EAe A )/D 

A 3 7 

= 

EI y , 0'(EAI<^ e A -EBJ)sin0/D 

A 38 

= 

-A yj cot 0 

A 39 

= 

EI y ,0'(EB* e A -EI Z , k 2 A )/D 

A 43 

= 

-G3 e A T q (EI Z , -EA e*) sin 0/D 

A 46 

= 

A 38 

A 47 

= 

G3 (EI Z , -EI yt -EA e A ) sin 0 cos £ 

A 48 

= 

~A^ 3 sin 0 /(e A T q ) + G3 EI y , cos' 

A 49 

= 

G3 EI yl e A cos0 /D 

A 53 

= 

-A^ 3 cot 0 

A 56 

= 

-a 37 

A 57 

s 

A^ 3 cos 0 cot 0/(e A T q ) -G3 EI y 

A 58 

= 

“ A 47 

A 59 

s 

A 49 tan 0 
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A 61 s 

me (n Z + to z ) sin0 

A 62 s 

2 

-me a) cos 0 

A 63 s 

-mk 2 w 2 + mfl 2 (kf 1 -k 2 ) cos 
m iii 2 * • 1 1 

A 6* " 

2 

-m n ex sin 0 -(e^ T Q sin 0 )' 

A 65 = 

m ft ex cos 0 + (e^ T Q cos 0 )' 

A 7I s 

-2 im we sin 0 

A 73 = 

2 

-m n ex cos e 

A 75 = 

- T o 

A 7,U s 

1 

A 81 = 

-A^| cot 0 

A 83 = 

n 

-m ft 4 ex sin 0 

n 

■a- 

OO 

< 

T o 

A 8, 10 = 

-1 

A 91 * 

-2 im flw 

A 92 = 

-m a) 2 3 

pc 

A 10,l = 

-m (fi 2 +03 2 ) 

A 10,2 s 

-2 im ftw(3p C 

A 1 0,3 = 

A 61 

A 10,4 = 

~ A 81 

> 

jp 

ii 

A 71 

a j.i,i = 

" A 10,2 


€ 


f 
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where 


**. *.t*h K*\ 


A ll,2 = 

-mu 

A ll,3 s 

A 62 


D = G3EI y ,(EI z ,-EAe^) 


PENDULUM EQUATIONS 

The pendulum configurations considered are of the simple type which are free to 
oscillate in the flap and lead-lag modes. The free vibration equations of motion for 
simple pendulums in both the uap and lead-lag degrees of freedom are developed in an 
identical fashion. These equations are formulated by first determining the acceleration 
of the pendulum mass in the inertial system and then transforming it to the deformed 
blade system. The absolute acceleration is used to express the inertial reaction force 
and moment at the pendulum hinge. Setting the reaction moment about the pendulum 
hinge to zero yields the equation of motion. These nonlinear equations of motion are 
subsequently linearized by assuming small oscillations for the pendulum about the 
steady-state condition. By assuming a simple harmonic solution for the linearized 
equation of motion of the free pendulum, the deflection of the pendulum is obtained in 
terms of the blade deflection. The point transfer matrix at the pendulum hinge station 
on the rotor blade is derived by simple equilibrium considerations of forces and 
moments across the pendulum after linearizing the inertial forces and moments. 

Acceleration Components 

The acceleration of the pendulum mass is determined with respect to an inertial 
reference system. This entails appropriate transformations from the inertial to the 
rotating system, to the preconed system, and finally to the deformed blade section. The 
derivation of the acceleration components is given in Appendix B. Equations (B-18) to 
(B-20) for these acceleration components are reproduced here as 


**•’!? • » 
a x , S 5 ~2G v - n (x + 5 ) + 2 n c Sine -2 £1 q cos 0 (24) 

a , = ri + (Cl x v + w + ft 3 x) sln6 +(v+2ftE; -O^v 
y pc 

+ fi 2 x u) cos 0 - & 2 n cos 2 e + ^ ft 2 C sin 2 0 (25) 
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i 


a , s ? -(V + 2 ft? - n 2 v + fi 2 x p) sin 0 + ( fl 2 xV + w 

M 


+ fl 2 3 pc x)cose - fl 2 £ sin 2 0 + ^ n 2 n sin 20 (26) 


where (£,r),0 is a general point relative to the deformed blade coordinate system. It 
may be noted that the gravitational contribution to these accelerations has been 
neglected. 


Flapping Pendulum 

The blade-mounted pendulum absorber is shown In figure 3. The pendulum after 
deformation is shown in figure 4. The deformed coordinates of the pendulum hinge and 
mass can be related from figure ft and are given by 

E = A cos A 
P 

^p = n A 

5p = C A + A sin A 

where 

A 5 the pendulum arm length 
A 5 rotation about (0, ti a , C A ) In the x' z 1 - plane 

| = 0(c) ; A =0(1) 

Differentiating equations (27) with respect to time yields 

tp = “A A sin A (28a) 

V* tt • o 

E = -A A sin A - A A cos A (28b) 

r 
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n\ 


% 

5 0 

( 28 c) 

"p 

S 0 

( 2 8 d) 

'p 

• 

s Z A cos A 

( 2 8 e ) 

M 

«» 

a fl( A cos A - A 2 sin A) 

( 2 8f ) 


Substituting the above expressions into equations (2^) to ( 26 ), the acceleration 
components to the first order (second and higher order terms are discarded) for the 
flapping pendulum become 

a^, = - % A sin A A 2 cos A - 2 ft v A - ft 2 x A 

- ft 2 % cos A + 2 ft S, A cos A sin 0 ^ ( 29 ) 

a y , • ( a 2 x A V A + w A + R 2 $ PC x A ) sin 0 A 

, 2 2 

+ (v A -2 ft & Asin A ft v A + ft x A y A ) cos 0 A 

- ft 2 q A cos 2 0 A + ^ fi 2 ( 5 A + Jl sin A) sin 2 0 A 

( 30 ) 

** * 2 ♦ 
a , = fi, A cos A - l A sin A 

2 p 

* 2 2 

- (v A -2 ft £, A sin A - ft v A + ft x A y A ) sin 0 A 

+ <si 2 x A v A + w A * fl 2 p pc * A ) c ° s e A 

2 2 12 

- ft ( C A + Jl sin A) sin 9 A + ^ n A sin 2 0 A 

( 31 ) 
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These acceleration components will be used to express the inertial reaction force and 
moment at the pendulum hinge. t 

The total inertial reaction force acting on the pendulum mass, M, can be written in 
terms of the above acceleration components whose associated unit vectors are i', jYk' as 





a,, k') 


( 32 ) € 


The position vector of the pendulum mass with respect to the hinge point is 


t = Z (cos A*f' + sin aIc') 



( 33 ) 


and the inertial moment reaction about the hinge point is 


or 


M = *r x £ 


r 


k' 


= -M & 


cos A 


sin A 


V 


V a 2' 

P 


M 


Mia , sin Ai' 
y P 


+MA (a , cosA-a , sin A) j* 

2 X 

P p 


-M 11 a , cos A k' 
Y P 


( 34 ) 


In addition to the inertiai reaction moment, there is a moment due to damping at 
the pendulum hinge. This damping will be presumed to be of the linear viscous type. 
The corresponding moment reaction about the hinge can be represented as 


20 


fii A = cnA|' 


(35) 


where 


N 


= steady-state normal force on the bearing due to F Ax , 


and F Ay , 


C 5 pendulum hinge friction coefficient 


Since the flapping pendulum is free to rotate about the y'-axis, the y' component of 
the inertial moment reaction is in equilibrium with the damping moment. Using 
equations (34) and (35), the equilibrium equation can be written as 


MH(a x , 


sin 


4-a . 


cos a) -CN A = 0 


(36) 


Equation (36) provides the pendulum equation of motion. The other two components of 
the inertial moment, the damping moment, and all three components of the force 
reaction provide the pendulum reactions which act on the blade at the hinge point as 


F Ax' 

= -Ma , 
X P 

(37) 

F Ay’ 

= -Ma , 
y P 

(38) 

F Az’ 

- -Ma z , 

(39) 

M Ax’ 

= MU a , sin A 
y P 

(40) 

M Ay' 

s -CN A 

(41) 

M Az- 

= -MHa , cos A 

Y n 

(42) 
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Substituting for the acceleration components in equation (36), the pendulum 
equation can be written as 

#6 9 9 

% A' - v A + SI x A P A ) sin0 A cos A 

+ (n x A v A + w A t n b p x A ) cos 0 A cos A 

*c 

2 2 

+ Q & cos Q a cos A sin A 
+ (2flv A + fi 2 x A )sinA 

+ (- fi 2 5 a sin 2 0 A + j fi 2 n A sin 2 0 A ) cos A s -CN A/M& {it 3) 

It may be noted that in the above equation the rotary inertia of the pendulum is 
neglected. This implies that the pendulum consits of mass, M, which has very small 
dimensions (point mass) in comparison with the pendulum arm. 

The nonlinear equation of motion given by equation (43) is linearized for smali 
angular oscillations, 6> about the steady-state equilibrium angle of the pendulum, 
where 

A = <S + 6 0 (44) 

Since the perturbation angle, 6 , will be fairly small as 

$ 2 « 1 

the following approximations can be used. 

cos A = cos 6 0 - 6 sin <S (45) 




i 


if. 


sin <$ s sin <$ 0 + 6 cos 6 Q (4 6) 

The steady state equilibrium angle, 6^* is obtained from equation (43) by setting 
the acceleration and velocity terms to zero°(assuming that the elastic deformations are 
zero). This process yields 

X A tan ^o + cos2 0 A^ sin 6 o s ^A sin2 0 A “1 h A sin 2 0 A - 3 x A cos 0 A (47) 


In the above equation 0 A is the steady-state value of the pitch angle at x A . It may be 
noted that the steady equilibrium angle, <$ Q , is independent of the rotor speed. This 
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result is due to neglecting the gravitational force in equation (36). The effect of the 
gravitational force on the steady equilibrium angle and the uncoupled pendulum 
frequency is very small (ref. 7). 

Substituting equations (45) and (46) into equation (43) and subtracting the steady- 
state relation, equation (47), from the resulting equation, the following equation of 
motion for the perturbation pendulum angle is obtained. 


U 


+ CN t 
+ M S, 6 + 


ir 


6 ( : 


cos 



slns 0 


+ 



V A + w A ) cos e A 


<V°‘ 


V A + 


n 2 X 


A 


y A )sln 0 a] cos 6 o 


+ ficos 


S o- 


[<« 


A 


V A + 


w A ) cos e A 


-(v A - fi 2 v A + fi 2 x A P A ) sin e A ] sin 6 Q = 0 


(48) 


In equation (48) the 6th and 7th terms (underlined terms) are nonlinear. Comparison of 
the 7th term with the 5th term illustrates that this nonlinear term is of order 6 tan o 0 
compared to unity for the linear term. Since both 6 and <$_ will be fairly small then 
the 7th term can be rationally neglected. Comparison of the nonlinear 6th term with 
the linear 4th term indicates that they are of the same order of magnitude since 



Both of. these terms will be fairly small. In the interest of linearization and the desire 
to retain an explicit damping term, the 6th term will be dropped and the 4th term 
retained. 

The steady-state normal force on the bearing can be written as 

N ■ OL' +P Az '] 1/2 

where F A „ and F. 
may be noted thatTT 

the steady-state value. Substituting for the acceleration components (eqs. (29) and (31)) 
in equations (37) and (39) and retaining only the steady-state terms, equation (49) 
becomes 


„ are the steady-state values of F A „ and F A „ respectively. It 
is generally a time-dependent force' but is approximated here by 


N = 


2 

M £2 (x A ^ 


ft cos 


«o> 


(50) 
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It may be noted that high order terms are neglected. Substituting equation 00) Into 
equation (48), the linear pendulum equation of motion becomes 

I 

X X \ 

& 6 +cn 2 (^ + cos 6 0 ) 6 + + a cos 2 0 A cos 2 6jS 

€ 

+ 2flv A sin 6 Q + (w A + ft 2 x A v A ) cosO A cos6 Q 

“ ( ^A " q2 v A + 0,2 X A sIn °A cos 6 0 c 0 (51) 

The uncoupled pendulum equation can be obtained from equation (51) by dropping 
the blade deflection terms. The resulting equation can be written as 

6 + 2&0) p 6 + w 2 6 s 0 (52) 

where 


“p ' n («^T57 + cos2 s o cos2 e ft) 1/2 


Cfl(x A /fc +cos6 0 ) 


(53) 

( 54 ) 


It may be noted that the pendulum hinge coefficient, C, in equation (54) can not 
normally be specified, but representative values of the damping ratio, g , can be 
estimated. 

The inertial reaction loads applied to the blade due to the pendulum will be 
considered to act at the elastic axis pendulum hinge location, (x A ,0,0). These reactions 
will be treated as discontinuities in the tension, shear, and bending moment distributions 
along the blade. The values of these discontinuities can be written in terms of the 
component reactions as 


c* 


^ St Vi#l 


* 


1M x- 

— 

M Ax' + 

r, A F Az» - C A F Ay' 

(55) 

‘M yl 

S 

M Ay ,+ 

C A F Ax' 

(56) 

AM Z , 

= 

M Az' " 

n A F Ax' 

(57) 

AV x' 

= 

Pax- 


(58) 

A V 

a 

Pa/ 


(59) 

AV z' 

= 

F Az' 


(60) 


Substitution of equations (29) to (31) and equations (37) to (42) into the above equations 
yields 

AM x ,/M = (ft sin A + S A ) [(^ 2 X A v a + + n 2 ^pc x A^ sin e A 

+ (v A ~2 ft ft AsinA-ftS^ + Jrx A y A )cos0 A 

2 2 12 *1 

- ft n A cos 6 A + ^ ft ( C A + ft sin A) sin 2 0 A J 

- r) A [ft A cos A- ft A 2 sin A-(v A ~2ft ft A sin A 

2 2 2 

-ft v A + ft x A y A ) sin 0 A + (ft X A v A + w A 

+ ^2 3 pc X A ) cos 0 A - ft £ A sin 0 A 

- ft 2 ftsin Asin 2 e A + 5 n 2 h A sin 2 °A 3 ( 6 ! ) 
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A M ,/M s ^ [a A sin A + i A 2 cos A+ 2ft v A + ft 2 x A 
+ ft 2 & cos A ~2ft A A cos Asin 0 A ] 

-C ft 2 (x A + a cos A) A 

A M z ,/M s ~ n A A sin A + it A 2 cos A + 2 ftv A 

2 2 • "i 

+ ft x A + ft £cos A -2 ft H AcosAsin 0 A J 
r* 2 d, 2 

- jlcos A [_( ft X A V A + + ^pe x A^ s * n 0 A 

• O o 

+ (v A **2 ft % Asin A- ft v A + ft x A p A ) cos 0 A 

- ft 2 n A cos2 0 a + i ^ 2 (t a + ^ s * n a ) sln 2 e A ] 

A V^,/M = Jl A sin A + % A 2 cos A+ 2ft v A 

+ ft 2 x A+ ft 2 & cos A-2 ft H A cos Asin © A 
A V y ,/M = -( SI 2 x ft v A + w A + a 2 3 p0 x A ) sin 0 A 

- (v A -2 ft A Asin A - ft 2 v A + ft 2 x A y A ) cos © A 

2 2 12' 

+ ft " n A cos 0 A ft ( C a + ^ sin A) sin 2 0 A 


( 62 ) 


( 63 ) 


( 64 ) 


( 65 ) 


( 66 ) 


A V z ,/M = - % Acos A+ i A 2 sin A + (v' A -2 ft i A sin A 

2 2 2 
- ft v A + ft x A n A ) sin 0 A - ( ft x A v A + w A 

+ ft 2 6 pc x A ) cos0 A + ft 2 2 a sin 2 0 A 

+ ft 2 il sin Asin 2 0 A ft 2 n A sin 2 0 A 

where M is the pendulum mass. 

Using the adopted ordering scheme given by Table 1 and equation (28), the orders of 


magnitude of the above expressions are 

AM X ,/M 

as follow 

fW 

0(e 2 ) ' 


AM yI /M 

T** 

0(e) 

> (67) 

AM ,/M 
z 

IV 

0(e) , 

1 

AV X ,/M 

f>» 

0(1) \ 

1 

4V y ,/M 

I’M 

0(e) | 

> 

> (68) 

AV Z ,/M 

A* 

0(e) ) 



For inclusion in the blade equations, the shear and moment reactions must be 
transformed into the undeformed x, y, z system using the transformation matrix of 
Appendix B. The transformations consistent with the above orders of magnitude are 


A M x /M = ( AM x ,/M) - (A M y( /M) ( p A cos 0 A + v A sin0 A ) 


+ ( AM z ,/M) ( p A sin 0 A - v A cos 0 A ) (69) 

A M y /M s ( A My,/M) (cos 0 A - ^ sin 0 A ) - ( A M z ,/M) (sin 0 A + <J> A cos 0 A ) (70) 

AM z /M = ( AMy,/M) (sin 0 A + $ A cos e A ) + ( AM z ,/M) (cos 0 A - <|> A sin 0 A ) (7 1 ) 
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2 

where only 0(e) terms and larger are retained. Similarly 

« 

AV/M » A V ,/M (7 2 ) 

X X 

A Vy/M S (AY X ,/M) y A + (AVy,/M)eos 0 A -(AV 2 ,/M) sin 0 A (7 3 ) 

A V z /M s ( A V x ,/M) V A + ( A V ,/M) sin 0 A + ( A Y^/M) cos 0 A ( 7 ‘I ) 

where only 0(e) terms and larger are retained. 

Substituting equations ( 61 ) to (66) into equations ( 69 ) to ( 74 ) and retaining only the 
linear time dependent terms (to be compatible with the completed time dependent 
blade equations), the shear and moment reactions in the undeformed preconed system 
are obtained after making use of equations ( 44 ) to ( 46 ). 

A M x /M = [ U sin + ? A ) sin 0 A - n A cos 0 A ] w A 

+ [(£ sin c a ) cos 0 A + n A sin0 A 6 sin 6 q 

2 2 2 1 
-ft v A + ft x A y A + ft £ <S cos 6 q sin 0 A J 

+ 9, 6 cos <S Q [ ft 2 3 p C x a sin 0 a ~ ft ^ n a cos 2 0 A 

+ ft 2 ( + i sin fi 0 ) sin 0 A cos 0 A ] - n A & 6 cos 6 q 

- n 2 x A ( n A sin 0 A + ? A cos 0 a ) y A 

+ S, sin 6 q ft 2 x A v A sin 0 A (75) 
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A My/M S (n A sin 0 A + C A cos 0 A ) [a 6 sin S Q + 2ft v A 
-ft 2 A 6 sin 6 q -2 ft A 5 cos 6 q sin © A ] 

+ Acos6 0 sin 0 A [( fi 2 x A V A + w A ) sin © A 
+ (v A - 2Q& <5sin S Q -* S^S 2 v A + ft 2 x A y A ) cos 0 A 
+ ft 2 A 6 cos 6 q sin 0 A cos 0 A J -ft 2 £ 6 sin S Q sin 0 A 

• Opc X A sin 9 A - \ cos2 9 A + < «A 

+ A sin 6 q ) sin 0 A cos 0 A j + ft 2 x A ( ri A cos 0 A 

- C A sin 0 A ) <j> A -C ft 2 6 (x A + A cos <5 q ) cos 0^ 

A M z /M = (- n A cos 0 A + C A sin © A ) [a 6 sin S Q 

+ 2 ftv & ~ ft 2 A 6 sin 5 -2 ft A 6 cos 6 sin 0- 1 

. A o o A J 

- A cos SqCos © a [( ft 2 x A \> A + w A ) sin © A 

+ (v A -2ft A 6 sin 6 Q - ft 2 v A + ft 2 x A y A ) cos 0 A 

2 "i 

+ ft A o cos <S Q sin 0 A cos 0^ J 
+ ft 2 A 6 sin 6 q cos e A [ 9 pc X A sin 9 A 

- n A cos 2 0 A + ( S A + & sin <$ 0 ) sin 0 A cos 0 A ] 

* 

V ■ 

2 

+ n X A ( n A sin 0 A + C A COS 0 A ) <|> A 
-C ii 2 <S(x A + A cos 6 q ) sin 0 A 


( 76 ) 


( 77 ) 
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A V x /M s sin 6 q + 2 fl v A - n 2 l 6 sin S Q - 2 to l <5eos S Q sin 0 A ( 7 8 ) 

A Vy/M - -v A + 2 ft *, 5 sin S Q + fl 2 v A 


4 


- n' 


6 cos 6 q sin 


0 A + 


H 6 cos 


& sin 0 A 

O A 


(7 9) 


A V z /M a -w A - A 6 cos $ 0 cos 0 A (80) 


Since these discontinuous ioad reactions contain the time dependent pendulum 
angle, it must be related to the elastic displacements for its elimination. This can be 
accomplished by using the linearized pendulum equation (eq. (51) . By seeking the 
simple harmonic motion of the form 


6 (t) = 6 e ia) t 


(81) 


and corresponding harmonic solutions for the state variables, the pendulum deflection 
can be written as 


where 


6 ‘ v A 5 v + w A 6 w + + Vv 


S v = - [( U 2 + ui 2 ) cos sin e A +2intusin 6 Q 3 /a 

= w cos <5 cos 0./a 

W O A 

x. cos 6 sin 0 A /a 
U A o A 

9 

5 v 58 " ^ X A cos o C0S 9 A^ a 


( 82 ) 


( 83 ) 

( 84 ) 

( 85 ) 

( 86 ) 
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a = & (w 2 - w 2 ) + 1C fi 2 a) (“— + cos 6 q ) 


(87) 


where C can be determined from the estimated damping ratio of equation (54) as 


C a 


2 *e<idrr + cos2 s o cos2 6 a )1/2 

o 


S2 (x A /jl + cos 6 q ) 


The discontinuous load reactions associated with the presence of the pendulum are 
used to relate the bending moment and shear force state variables on the right side 
(outboard) of x A to those on the left side (inboard)* These state variables are indicated 
in figure 5. The state variables are related as 


M 

M 

M 


M L - AM 

v R - 

V L 

- AV 

X X 

X 

X 


m l - am 

V R = 

V L 

- AV, 

y y 

y 

y 


M^ -AM 
z z 

V R - 
z 

v L 

z 

- AV 


( 88 ) 


I 

) 


All other state variables must be continuous across the point x«. 

Substituting for the shear and moment reactions, given'by equations (75)to (80), 
with the explicit time dependency, the above relations can be described by 


: z R j = [p] jz L J (89) 

where [P] is a "point" transfer matrix. The elements of [P] are all zero except for 


P.. = 

ij 

1 

(i = j) 


P 61 ' 

M £($, sin 6 q + 

C A ) cos e A + n A sin 

e A ] (ft 2 + w 2 ) - 6 y M 

P 62 = 

M [( i sin + 

C A ) sin e A - n A cos 

6 a]“ 2 - { w M x6 

P 64 ■ 

-M H sin 6^, cos 
0 

1 9 A fl2x A-% M X 6 


P 65 ■ 

-Mfcsin 6^ sin 
0 

0 A n 2 x A -6 v M x6 



x 6 
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-M [( n A sin 0 A + C A cos © A ) 2iw ft 

-A cos « o sin0 A cos 0 A (w 2 + ft 2 ) ] - <$ v M yg 

Mw 2 A cos 6 Q sin 2 e A - M y6 

2 

-M ( T 1 A cos e A - c A sin 0 A ) ft x A 
2 

-M ft x A % cos S Q sin 0 A cos 0 A - <$ ^ My 

-MiJ 2 x A «.co S 6 o sin 2 0 A -6 v M y{ 

-M cos C Q cos' 2 © A ( to 2 + ft 2 ) 

-( n A cos 0 A - ? A sin 0 A ) 2i a) ft J - S v M 2 6 
2 

-Mto A cos <S Q cos 0 A sin © A ~ <S W M z ^ 

2 

-M ( n A sin 0 A + C A cos 0 A ) ft x A 

M ft 2 x. A cos cos 2 0. - 

A o A y zo 

M ft 2 x A % cos <$ 0 cos 0 A sin 0 A - 6^ M z ^ 

-2 M its ft - 6 V . 

v x 6 



$ 


P 95 s 

■ 5 v V x 6 

I! 

o 

a. 

-m(oj 2 + n 2 )- s v v y{ 

P 10,2 5s 

- S w V y 5 

p 

10,4" 

-% v y6 

P 10,5 “ 

- 5 v V y 6 

p ii,i“ 

' 6 v V z S 

p ii,2 s 

- m “ 2 - 5 w V z« 

P ll,4 = 

- V*« 

P ll,5 = 

- 6 v V z « 


where 


M x 6 = M [ (jUn6 o + ?A )cos 9 A + n A sin 0 A] (fi cos 5 o sin 0 A 


-ltd sin 6 ) 2 ft &- Jl n A cos <5 ( ft^-w 

O AO 


,«r% 

+ ft 2 0 pc il x A cos 5 0 sin 0 A J (90) 


<rN 


33 


ii 



-M j ( r) A sin 0 A + cos 0 A ) it sin <S Q (w 2 + ft 2 ) 

+ [(A sin 6 0 + c A ) sin 0 A - n A cos ej 
2 

• ft A sin sin 0 A cos 0 A + (x A 0 _ sin 6 
o A A A pc o 

- itcos 2 6 q cos 0^) ft 2 itsin 2 0 A + [(it sin <$ 0 

+ C A ) cos 0 A + n A sin 0 A ] 2iu>ftitcos <S Q sin 0 A 
+ iC ft 2 w (x A + it cos 6 q ) cos 0^ | 

M j ( n A cos 0 A - c A sin 0 A ) Z sin <S Q (w 2 + ft 2 ) 

+ [(it sin 6 q + C A ) sin 0 A - n A cos 0 A j 

• ft 2 Z sin S Q cos 2 0 a + (x A 0 p C sin S Q 

-it cos 2 <S Q cos © A ) ft 2 it cos ® A sin © A + [( r)^ cos 0 

- ? A sin © A ) sin © A + Z sin <S Q cos 2 © a J 

• 2 i wft it cos S -iCft 2 w (x A + it cos 6 ) sin 0 A { 

o A o A ; 


(91) 


( 92 ) 




s -M [a, sin (to 2 + ft 2 ) + 2 i toft % cos 6 Q sin 0 A ] 

(93) 

v y 6 

5 -M [it cos S Q sin 0 A (to 2 + ft 2 ) -2 Ito ft % sin 6 Q ] 

(99) 

V, r 
z 6 

2 

= Mo) % cos 6 0 cos 0 A 

(95) 


Lead-Lag Pendulum 



The lead-lag pendulum after deformation is shown In figure 6. The equations of 
motion for the lead-lag pendulum can be derived by following the procedure used in the 
flapping pendulum case. This process yields the following equations for the lead-lag 
pendulums 

The relation between the pendulum mass and hinge in deformed coordinates 


Sp - & cos A 1 

rip = n A + i sin A \ 

? P ? A } 


Pendulum angle 


(96) 


A = 6 


(97) 


where 

6 q 2 steady-state equilibrium angle 
6 2 time dependent perturbation angle 
The equation for determination of the steady-state equilibrium angle, becomes 


x A tan 6 Q + % sin 2 0 A sin S Q = - A sin 2 0 A + ri A cos 2 0 A - 0 pc x A sin 0 A (98) 
The pendulum deflection in terms of the blade deflections 


6 


V A V 


w 


6 W + 




6 + 


v, 


(99) 
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where 


<S 

6 

<5 

6 


v 

w 

P 

V 


a 


w 


£ 


P 


C 


6 


[(ft 2 * to 2 ) cos 6 q cos 0 A -2itofl sin 6 Q ]/a 


n 

to cos 6 sin 0 A /a 


-f* x A cos 6 Q cos 0 A /a 


2 

- ft x A cos 6 Q sin 0 A /a 


& (to 2 -to 2 ) + i Cft 2 to (-— + cos 6 q ) 


- ft 


<3./ Xa rt 5 \ 

— -fr - r - + cos 6 sin 0 « ) 
cos 5 o A/ 


n/ 

Te 


2 “S 2 V 1 "" G a) 1/2 

~~ ft (x A /a + cos 6 q ) 


6 e i(ot 


The nonzero eiements of the point transfer matrix become 
p i) • 1 

P 61 s M [^ A cos 0 A + ( + & sin S Q ) sin 0 A ] (to 2 + ft 2 ) - 

P 62 = M A sin e A - ( 1 A + * sin «o) cos 0 A ] 0) 2 - « W M X g 

P 64 = " M S o 5in 9 'A 0 2 x A - 6 p M X 5 


( 100 ) 

( 101 ) 

( 102 ) 

(103) 

(104) 

(105) 

(106) 
(107) 


36 


P 6J = Mjlsln « 0 cos 0 A f! 2 x A - « V M X{ 

P 7i s -M [( n A sin 0 A + ? A cos G a ) 2iwft + % cos <5 Q cos0 A 

♦ sin 0 A (w^+ ft^)[| "^v^y6 
P 72 = MU cos 6 q cos 2 0 A w 2 - 6 W M y fi 

p 73 = -m ( n A cos e A - c A sin e A ) a 2 x A 

P 7 d « Milcos 6 q cos 0 A sin 0 A Q 2 x A - 6 y M y 6 
P 75 = ~M A cos S Q cos 2 0 A ft 2 x A - <S V M y 6 
P gl = -M [$, cos 5 o sin 2 6 A ( 0 ) 2 + ft 2 ) 

-(n A c °s 0 A - S A sin 0 A )2iwft] - <$ y M z 6 
P g 2 ~ M to 2 il cos <$ 0 sin 0 A cos © A “ $ w M z $ 

P g3 = -M ( r\ A sin 0 A + C A cos 0 A ) ft 2 x a 

P 8if = Mft 2 x A Acos6 o sin 2 0 A -6 y M z6 
P g5 = -M ft 2 x A ft cos 5 0 sin 0 A cos © A - <$ v M z 
Pc^ = -2Mi toft - <5 y V x ^ 


where 


t> a» mm 


1 

10,1“ 

-M(u) 2 + n 2 ) 

10,2 s 

- 6 V 
°w v y 6 

’io,/r 

' 6 y v y 6 

' 10 , 5 s 

" Vy « 

11 , 1 s 

- { v V 2 6 

11,2 s 

-Mid 2 

11,4 s 

‘ 6 u V z6 

11,5 s 

- V* 0 


6 V * 
v y 6 



M x { = -M { [? A cos e ft + ( n A +«.sln « 0 )sin 0 ft ] 

* (8 cos 6 q cos 0 a + iw sin 6 Q ) 2 ft & 

-JU A cosS 0 (n 2 - w 2 )+ fl 2 3 pc X,x A cos 6 q cos 0 A J ( X 0 8 ) 

M y 5 s " M { ^ sin °A + r 'A cos 0 A^ ^ sin ' W 2 + 

+ [- S A sin 0 A + ( n A + ^ sin S Q ) cos 0 A ] tt 2 Z sin S Q sin 0 A cos 0 A 

+ (x A 3 pc sin 6 q - Zcos 2 5 q sin 0 A ) Q 2 Z cos 2 0 A 

“B n A + &sin <S o )sin0 A + cos 0 a] 2l wn Acos 6 o cos0 a 
- iCn 2 o) (x A + % cos S Q ) sin 0 A J (109) 
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M z 6 s M j ^ n A cos 0 A “ r *h sIn 0 A^ A sin 

+ [ C a sin 6 a "( n A + ^ sin 6 Q ) cos 0 A ] 8 2 1 sin S Q sin 2 0 A 

2 2 

-( x A 0 p C sln 5 0 cos fi 0 sin 9 a^ ^ &cos0 A sine A 
- [( n A cos 0 A - C A sin 0 A ) cos e A - Hsin S Q sin 2 0 A ] 

•2iwft % cos 6 0 -iC fi 2 w(x A + % cos 6 Q ) cos 0 A J (HO) 

V x 6 = “ M sin 5 o ( w2+ ^ " 2ia) Z cos 5 o cos 9 A ] Oil) 

V y 6 = M cos fi o cos 9 a ^ 2 + fi 2 ) + 2i ^ n A sin 6 0 ] (112) 

V z 6 s M w ^ cos *0 sin 0 A (113) 


AERODYNAMIC FORCES 

This section is concerned with the calculation of the aerodynamic loads on a 
helicopter rotor blade in forward flight. Due to the complexity of the flow field 
involved and the interaction of different effects in the flow, a completely general 
solution to the problem has not been obtained. There are a large number of 
approximate methods available. These methods range in complexity from simple blade 
element representations to lifting surface models with freely distorted vortex wakes, 
with associated ranges in computational expense, accuracy and detail of the solution. 
Commonly used is an approach based on a combination of simple momentum theory and 
the blade element description. This method is adopted for the present study. The 
method takes no note of the reverse flow and stall. 

The aerodynamic forces are formulated from two-dimensional, incompressible, 
quasi-steady, strip theory in which only the velocity components perpendicular to the 
spanwise axis of the deformed blade are assumed to influence the aerodynamic loading. 
Account is taken of the varying freestream velocity associated with the rotating blade 
by employing Greenberg's extension (ref. 11) of Theodorsen's unsteady theory for 
determining the aerodynamic lift and pitching moment acting on the blade. The 
resulting expressions are specialized to the case of quasi-steady flow by setting 
Theodorsen's circulation function to unity. Simple momentum theory is used to 
calculate the steady flow induced by the rotor. The expressions for the aerodynamic 
loading are based on the analysis reported in reference 12. 
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In Theodorsen's theory a two-dimensional airfoil is assumed to be pivoted about an 
axis which may be distinct, in general, from the aerodynamic center. The airfoil is 
pitched at an angle e (t) to the freestream flowing at constant velocity V. The airfoil is 
vertically displaced with velocity h (t) as shown in figure 7, Greenberg (ref. 11) has 
extended Theodorsen's theory for a varying free-st.ream velocity V(t). The lift and 
moment acting on an element section of the blade may be expressed in terms of the 
circulatory and noncirculatory components as 


L = 

L c + l nc I 

in*) 

M = 

M c + m nc ) 


With the airfoil pivot axis (analogous to the rotor blade elastic axis) at the airfoil 
quarter chord (the airfoil aerodynamic center) the relations for the aerodynamic loads 

per unit length may be written as (ref. 

12). 


l nc 

= 2 p a T *“'^p + ? p * 

(li-5a) 

L C 

= i pacu(-u p + | e) 

(1 15b) 

m nc 

= ac (|) 2 (-Up -Uj e + ) 

(115c) 

M c 

= -|pac(f) 2 2U T e 

( 1 1 5d) 


The profile drag force acts parellel to U and is given by 

D = i p ac U 2 (U6 ) 


The components of the aerodynamic forces in the deformed coordinate system can 
be written as 


where 


F , = -L_sin a -D cos a 

y ^ 

F 2 , = L c cos a + L n c -D sin a 


(117) 


sin a = U p /U 
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cos a - Uy/U 


U 


- \[ u 


2 + u 2 

T + U p 


Substituting equations (115a), (115b) and (116) into equations (117) yields 

F y ■ 5' )ac [ u p-f u p fe - f r u T] | 

F z ' ■ fpacj;-U p U Tt feU r fU p+ (|) 2 e] j 
From equations (115c) and (115d) the pitching moment can be written as 

M K , <=M ) = -± p ac {f [u T c-U e] 


where 


(118) 


(119) 


Up, U^. = relative velocity components in the deformed coordinate system (fig. 7) 
e = angular velocity of the blade about the x' - axis 


V(t) a 4 U T + U p 

The assumed orders of magnitude of the quantities used in deriving the aerodynamic 
loads are 

U 

-rp = 0(e) 

U T 

^ . o (e 2 ) 

a 

I . o(«) 

srk = ° (e) 

0 = O (e ) 



The aerodynamic forces in the undeformed coordinate system are obtained by using 
the transformation given by equation (A-4). To second order, the aerodynamic ioads can 
be written as 


l v <=i y 


Fyl - F Z I ( 6 ♦ 


L w <= F z> 


Fyj (0 + <j>) + F z r 


(120) 




M x’ 


In equations (118) and (119) the blade aerodynamic loads are expressed in terms of 
Uj, Up, and e . In order to use the expressions for blade aerodynamic loads in the 
blade equations of motion, Uy, Up, and e must be expressed in terms of the blade 
bending and torsion deflections v, w, and cj> . The blade velocity is easily expressed in 
the x, y, z coordinate system. The deformed blade coordinate transformation, which 
relates the x,y, z and x', y', z' coordinate systems, is given by equation (A-4). This is 
used to express the Uy and U p velocity components in terms of v, w, and <j> . 

The total flow velocity seen by a point on the elastic axis can be written from 
reference 12 as 


v = ( ji f n r cos \}j + n r x e pc -b + n v) i 

+ ( w -v “ft{x + u) - y, ft R sin \jj)t 

IT V-. 

+ (ftRA - y f ftR3 pc cos ~w - ^3 pc v) k (121) 

where 

V cos a 
_ s 


V sin a -Vj 
A ft R 

and the aerodynamic velocity components are shown in figure 8. 

The tangential and perpendicular velocity components, Uy and Up, are obtained 
from equation (121) using £t]] from equation (A-4). To second 1 order in the dependent 
variables they have the form 
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’ If 



ft 


ji f fiR v' cos i]j + Ox+ |ij ft R sin tp + v (122) 
y^ftR w 1 cosi|j- y f HR v'0 cos^ 

■( y f ft R siru{J + ftx) (0 + <{>) -v0 


+ y f ftRB cos\|) -ftRX + w + ftB v (123) 
z pc pc 

It may be noted that the nonlinear terms are neglected. 

The angular velocity about the x’ axis, e , is now considered. It can be regarded as 
being composed of three parts: the first part due to the rigid-body angular velocity of 
the hub in space, the second part due to the control inputs, and the third due to the 
angular velocity associated with the elastic deformations. The total sectional pitching 
velocity can be written as (ref. 12.) 


e = (3 pc + w' + 0j c sin tp - 0 Js cost) + t (12 4) 


It may be noted that the total geometric pitch angle, 9 , is given by 


9 = 0 pt + 0 c -9 lc cos t -9 ls sin t (125) 


where 0 is the built-in twist angle (pretwist), 0 is the collective pitch angle 
measured Sf the tip, and 0, and 0 . are the first harmonic cyclic pitch components. 
The built-in twist is a linearlinear function of x and may be written as 





- e t (R - x) 


(126) 


The induced veioctiy, v-, is calculated by equating the integral thrust to the thrust 
from momentum theory using the relations 
S 


v. 

i 



(127) 


£ 
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1 


(128) 


' T ' 2 " ipli 2 ^ 


2,jt „ R 

// 


L w dxd iJj 


where L is the steady-state value of L . 

Substituting equations (118) and (119) into equations (120) and using expressions for 
Uy, Up, and e, given by equations (122) to (124), the aerodynamic loads to second order 
can be r written as 


L y = ^ p ac £-(ft x + pj ft R sin \|i)0 + 2 ( p f ftR0 p C cos \}> 


-ft R X)J i}/ + |pac £-(ft x + p f ftRsinty) 

• ( p^ ft R 3 p C cos i|>» ft R X ) J cp + - p ac j^-( ft x 
+ PjftR siniji ) 0 + 2 (p jft R 3p c cos \\> -ftRX )J 
• P^ ftR cos ip v + ^ p ac £-(ft x + p f ft R sin tp) 


>(p,ftR3 cos4> -S2 RX) 0 - — (ftx 
' f pc a 


+ p^ ftR sin ip )^ + ( ftR 3p C cos -ft RX )^ J (129) 
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j p ac £-(0 x + y^ 0 R sin ty) 0 0 pc J v 
+ ~ pac [2 (fix* y f ORsimjO 0 -(y f RR$ pc cos \p 
~ORX)] V + :| p ac £-(0x + y f ORsM)w-|w] 

1 r 2 

+ -pac[-(Ox + y f 0 Rsin 4 >) 

+ | y f O 2 R cos i|>] 4> + ^pac[^(fix 
+ y f 0 R sin 40 J j> + 5 pac £ 2 ( 0 x 
+ y f flRsini|>)e -(y f ORB pc cos \Jj -0 RX )] 


•(y f OR cosi|0 y + ^ pac -(Ox 


+ y^ 0 R sin \J j ) y^ OR cos 4> + | ft (0 x 
+ y f 0 R simp ) + | y f ft 2 R sin \p J v 
+ |pac y f 0 R cos^] v 

+ ^ pac [-( y f 0R 0 p(; cos^ -0 RX ) + ( y f 0 R sin ip 

+ 0 x) 0 ] ( Ox + y f 0 R sin \p ) 

+ |pac[|(O3 pc +|0)(Ox + y f OR sinty ) 

+ | y f 0 2 R cos 4)0 + f y f 0 2 R g pc sin \|i] 


( 130 ) 
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M ^ = - ~ P ac (|) 2 [-w + y f <2 2 R cos <i> 

+ 2 (ftx + ft R sin\|)) <j> + (ft 2 x + 2 y^ ft 2 R sin ij*) V 

- ja j ft R cos \>J 

- ^ p ac (|) 2 £ ( ft x + y^ ft R sin4> ) 

• ( ft p p C + 20 ) + y f ft R 0 cos vp 

+ y f ft 2 R 3 pc sin \Jj] (131) 

The above expressions can be rewritten in matrix form as (retaining only time- 
dependent terms) 

{fJ. [b (x, t)] | Z (x, t) | + jc<x,t)] jz (x, t)| + [d(x, t)] |z(x, t)} 

+ n Ej [ B c (x ’ l) ] n ( Z (x ’ x) } cos n 

"Ej [ C c (x - t) ]n{ Z(x ’ t) } cosn + 

3 

+ n E, [». (x > *>] n { Z (x ’ l) } sin n ^ 

+ ^ [C s (x, t)] n {z(x, t)| sinn 41 
+ n E^{ P c <x) } n cos n * + t, { P s (x) } n sin n * 


( 132 ) 



where j Z J is the state vector, defined by equation (1), and 

K} T = Lo o 0 0 0 -M 0 0 o-L v -L w J 


( 133 ) 


The nonzero eiements of the elastic-deformation-independent aerodynamic load 
vector, P, are: 

For n = 1} 


c6 

i p ac (|) 2 ft 

-2 ft x 6 | s + Uj ftR ( 0 t 

*'«>] 

n 

o 

u 

- i? pac ft R 

[-2f ipc |l f OJU -(0 t + 

6 c> 


+ iBpc“f 2 e Js -X« xe Xc ] 


'cl 1 = 

p ac ft R £ 

*-3 Pc y f ft x + f ( 0 t * 0 c) u f n 


2 

0 is"^ + i ^f 0 id 


s6 = 

|pac(|) 2 ft [ 

_ 2ftx0 ic + 2ft 0 pc p { r] 


slO = 

- ^ p ac ftR 

[-2~~ p f ftx+ (0 t + 0 c ) u f ft RX 


-Xftx0| s + ^B pc ii|ftR e lc ] 


si 1 = 

1 

- 2 p ac ft R 

[ ^ 0 p C ^U-f + Xfty|R 



2 

+ 2 ( 0 1 + e c ) u f ftx-(r|- +|u 2 ftR) e ls 


+ (r) e lc ] 
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For n = 2,* 


p c6 " 5 pac(f) 2 n Mj rr e lc ) 

r* C 

P clO = - j p acftR p 2 fiR + ^0p C ji 2 J2R 

+ ixv f nR 0j s + 5 s p0 ujJix e lc ] 

P cll = - j P ac nR [ "2 t + 0 C ) P j R R + R Pf x 0 i 5 

"I M f n 9 lc] 

p s6 “ |pac(f) 2 n (||i f aR e ls ) 

P s 10 ~ -? pacnR [-I (e t + 9 c> 9 pc^f fiR + ? 9 P c R f r ‘ x 9 ls 
-i X p f Si R 9 lc ] 

P sll = 4pacS!R[-ie pc p 2 nR-|cp f !! 0 IS -n M£X 0 lc ] 
For n = 3; 

P cl0 = - \ P ac (- | 3 pc n | U R 0 ls ) 

F^l | = ”2'P a c^F.(^ , jj|J2R 0 

P s io= -ipac nR(ie pc p 2 s!R e lc ) 

P sll= -^PacRRtip 2 S2R 6 ls ) 


The first seven terms in equation (132) represent the aerodynamic ioads due to 
elastic deformations. The coefficients of the corresponding matrices are given in 
Appendix C. 

The aerodynamic ioads as represented by equation (132) are to be included in the 
blade equations of motion. It may be recalled that the solution to these equations Is 
treated as individual contributions associated with each harmonic of the excitation 
loading. Such a solution may be written in a complex form as, equation (22), 

{ Z n (x,t)} = Re({z n (x)} e 1 **) (13») 

where 


(i) s nil (n = 1, 2, 3) (135) 

In the aerodynamic ioads expression (eq. (132)) the fourth to the seventh terms 
contain functions of the azimuth angle (rotational speed). Due to the existence of these 
terms, the above solution (eq. (134) can not be applied. Instead, the state vector should 
be written as 


1 7 {* 


t a_ 

V 


00 


D J V 

Yi=-o 


rn 




Jn t\ 


(1 


i c \ 
-'u/ 


In this case the resulting equations must be solved simultaneously for all harmonics. 
Even if only three harmonics are considered for the state vector (eq. (136)) obtaining 
the corresponding solution is very complicated. An approximate method is used to 
determine the aerodynamic forces by neglecting the fourth to the seventh terms in 
equation (132). These terms are of lower order than the corresponding terms in the 
blade equations of motion and thus may be neglected according to the adopted ordering 
scheme. Although the first three terms of equation (132) may also be neglected 
according to the ordering scheme, they are retained here. These aeroelastic terms will 
be the only representation of coupling between the elastic deformations and the 
aerodynamic loading. 

With the fourth to the seventh terms in equation (132) neglected the solution given 
by equation (134) can be applied. Substituting equation (134) into equation (132) yields 
the following complex notation for the nth harmonic of the aerodynamic load. 


{".} ■ Oa] H + M 


where 


{ ^a}“ am pii t ude vector of the aerodynamic load 

| p} = amplitude vector of the elastic-deformation-independent aerodynamic force 
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and 


{p} = 

{ p c}- 1 { p s } 

(138) 

K] = 

£b (x, ii) )J + iw jc(x,w)J - w 2 £d (x,ii) )J 

(13 9) 


From equation (139) and Appendix C, the nonzero elements of [A J are 

A a62 ' k p ac w 2 

A a63 = “ pac(|) 2 (2io)f2 x) 

a 1 _ /C\2 a 2 

A a65 = 2 pac y x 


A al0,2 = 
A al0,3 = 
A al0,5 = 
A al 1,1 = 

A all,2 = 
A all,3 = 
A all,4 = 


-^i pac w [-12x ( 0 c + 9 t ) + 2 Mf 6 ls " 2fiR * 

a 

- j p ac 12 2 x R A 

-JP ac(j fi 2 x jj f R 0 lc + fi 2 n 2 R 2 0 pc ) 

- ~ p ac (-12 2 x 0 pc ) 

- j ip ac 0 ) |^212 x ( 0 c + 0 1 ) ~ p f 12R 0j g + 12 R A J 

- j pac (-iio 12 x + | id 2 ) 

~~ pac(-12 2 x 2 -~ p 2 12 2 R 2 + q-p iiu 12 x) 
-ipac(-A M f R6 lc -i$2 2 M f 2 R 2 P pc ) 
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A all,5 s “ 5 pac(|ft 2 x) 


The induced velocity, v., Is obtained by substituting the steady-state value of L 
from equation (130) into equations (127) and (128). The resulting equation is 


where 


a n (6 X + C ) 
V 1 a H f tano> s -\ °-'y; T"7T 


fi^ * " 1 '' 


( 11 0 ) 


V cosa s 

= a r 


0 * 1 ) 


1 be 

a n 5 ^8 ffR a 


/ j /. n \ 


c n s (4 + 6 y 2 ) 0 C -(1 + 3 y 2 ) 0 1 R -6 y f 0 


Is 


+ |(33 pc +3 % e ic> 


(H3) 


Equation (140) is an algebraic equation in X • 


TRANSFER MATRIX METHOD 

o 

System Equations 

The blade equations of motion are given in a matrix form by equation (23) for the 
nth harmonic of frequency w as 

V-: 

{z(x)}’ = [aw] { Z(x)j + { r, (x) j ( 144 ) 

The boundary conditions corresponding to a hingeless rotor blade are idealized in this 
55 analysis as being cantilevered at the hub (x = 0) and free at the tip (x = R). 


51 


Hub; 

1 

z d (0)} B jo} 

(145) 

Tips 


(V R >i “ |°j 

(146) 

where 





1 g d <*>) T 

s Lv w ? li vj 

(147) 


M T 

= K M z V x V y vj 

(148) 


Equation (144) represents a set of li lingar first-order ordinary nonhomogeneous 
differential equations. It may be noted that [AJ is a complex matrix due to inclusion 
of the linear damping terms. The desired solution is obtained by the fomulation of a 
field transfer matrix which relates the state variables between any two points along the 
blade. 


Transfer Matrix 

The concept of the transfer matrix has been thoroughly discussed by Pestel and 
Leckie in reference 13. The backward transfer matrix [T(x) j is defined by 


Z (x)j = [t(x>] jz(o)} 


( 1 * 9 ) 


which relates the state vector at two stations. 

It is often possible to determine the transfer matrix by using simple statics. Such 
techniques are described in reference 13, These techniques prove satisfactory for 
simple lumped parameter or lower order uniform continuous systems, but become 
cumbersome when used for nonuniform and higher order systems because of the required 
algebraic manipulation. Murthy (ref. 9) has presented a systematic approach which 
eliminates much of the algebra and results directly in differential equations for the 
elements of the transfer matrix. The state vector usually satisfies the differential 
equation 


A | z (x) j = [a«] { z (x) 


( 130 ) 


Differentiating equation (149) with respect to x gives 
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dx 


z <*>) - i [H | z(0 >| 


(151) 


From equation (149) it is obvious that 


[ Z (0) j « [ t ( x )]' 1 jz(x) 


(152) 


Substituting this into equation (15 i) the following equation can be obtained 


ijz(x)| =|[tw][t«]-‘|z«) 


(153) 


Equating equations (150) and (153) gives 


F jk (v)l -A. Tt (v)1 

L dx L" K "'J 


r t 

L" 


U)~\ 


{z(x)j 


a ( 0 ! 


0 

( ) 


(15 4) 


Since equation (154) must be satisfied for ail values of x and all values of Z it follows 
that 


[ A <*>] ■ & [ T w] [ T <*>] ■' 

Then post multiplying both sides by [ T (x)] yields 

5 [ T «] ■ [ AW ] [ T(x> ] 


Therefore, the transfer matrix is given directly by the solution to equation (155). By 
letting x go to zero in equation (149) the initial condition of the transfer matrix 
becomes the identity matrix. 


[t<0>] - M “ 56) 

If equation (155) is solved as a coupled set of first-order differential equations, then 
equation (156) provides a sufficient number of initial conditions. 

The transfer matrix of the blade without a pendulum can be obtained by integration 
of equation (155) together with the initial conditions of equation (156). The tension 
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coefficient, T, appearing in matrix [A]is given by 


R 


T = 1 f mx^ dX| 


(157) 


When a pendulum is attached to the blade the transfer matrix at any spanwise 
location is obtained as follows. Let x A be the spanwise location of the pendulum hinge: 


i 

Case 1. 0 < x s x^ 


The transfer matrix in this case is obtained by the direct integration of equation 
(155) together with the initial conditions given by equation (156) similar to the case of 
the blade without the pendulum, but the tension coefficient is given by the following 
equation instead of the one given by equation (157). 


T 


n 


f mxj dxj + fVl i,'! 2 (x A + i cos 6 q ) (15 8) 


Case 2. 



< x < R 


Let [T.] be the transfer matrix i>f the blade at x s x^; this matrix can be obtained 
from Case 1 as described above. Recall that [[ P[] (see eq. (89)) is the point transfer 
matrix of the pendulum. Let [T 2 (x)] be the transfer matrix of the blade from x to x. 
This matrix is obtained by the integration of equation (155) with the initial conditions 
given by equation (156) from xf^ to x. While integrating this equation the tension 
coefficient given by equation (13/) should be used instead of the one given by equation 
(158). By the product rule of transfer matrices, the transfer matrix of the system at 
any spanwise location for this case is given by 


[t<x>] . [t 2 w] [p] [tJ 


Once the system transfer matrix is obtained, the state vector at any spanwise location 
can be determined. 


System Dynamic Response 


It has been shown how the transfer matrix can be used to solve the homogeneous 


differential equations. An application of such a solution is to find the natural 
frequencies and mode shapes of an elastic system (ref. 14). Once they are known it is 
possible to solve the most general cases of forced undamped vibrations, whether they 
are transient or steady-state in character, either for discrete or continuous systems, 
through the normal mode approach. On the other hand the steady-state response caused 
by a harmonic excitation is more readily solved with the aid of a particular integral of 
the nonhomogeneous differential equation without making use of the normal modes and 
natural frequencies (ref. (13)). 

The procedure adopted here for the steady-state response of forced vibrations with 
harmonic excitation is usually called the "exter ted" transfer matrix approach for 
undamped systems and the "complex extended" transfer matrix approach for damped 
systems. The procedure is valid for any linear damping since normal modes are not 
introduced. The restr ctions are that the response is steady and excitations are 
harmonic. With the knowledge of the system's response to harmonic inputs, the 
response to any arbitrary input can be computed by synthesizing the arbitrary forcing 
function from an aggregate of infinitesimal harmonic forces. 

The blade-pendulum equation of motion was given by equation (23) as 

jz(x)J ' = [a(x)] |z(x)J + { F a (x) | (160) 

By definition of the transfer matrix, the solution for the homogeneous part of equation 
(160) can be written as 


z h wj = [t (x)] {z<o) 


The particular solution of equation (160) can be written as 


M -/OH ft 


(s) 


ds 


(161) 


(162) 


where [T (x,s£j is the transfer matrix of the system from s to x (in contrast to [T (x)J, 
the transfer matrix from 0 to x). By the product rule of transfer matrices the following 
relation can be written 


[tw] ■ [ T <*■*>] [t«] 

From equations (162) and (163) the particular solution can be written as 

jfp(x)} = [TW] / [T (s )] - 1 { F(s)J ds 


(163) 


(164) 


The complete solution of equation (160) can be obtained by adding equations (161) and 
(164). 
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|z (x)J = [t (x) ] [z(0)j + [t(x|] / X [t(s)]‘ 1 {f(s)J ds (165) 

Equation (165) gives the amplitude vector of the simple harmonic response. The vector 
F % (0) is not completely known, but it can be determined from the boundary 
conditions of the system. ( _ , 

To determine the vector Z (0) , first rewrite equation (165) as 


where 


( z « <x) 1 


' T a <x> ; r b (xf 

( 

( Zj (X) ' 

. = 

_T C (X) j T d (x)_ 

1 


Z d (°) 


Zj (0) 



-T a (X) , T b W 

X 

r s a < s) | 

S b 

+ 

1 

f 

0 

■ + 


Jc (x) 1 T d (x)_ 

L S c (s) 1 

S d (s)j 


F d (s) 
> ds 


F f (s) 


(166) 


fz 1 T 
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0 
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[ T(s) ]"‘ 


The boundary conditions corresponding to a hingeless rotor blade are 

0 


z d ( 0 ) 


Z f (R) 


0 


(167) 

(168) 

(169) 

(170) 

(171) 


Substituting equations (167) and (170) into equation (166) yields 
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T b M 

T d W 



} 


'T a w ; T b w" 

r x 

X < s > 

_T c (x) | T d W 

J 

_s d « 



( 172 ) 


From equation (172) 



1 R 

"S b (sf 

+ [ 

T c (R > ! T d «>] f 
0 



{Ff(s)} 


ds 


(173) 


Substituting the boundary condition of equation (171) into equation (173) yields 
{ Z f (0) } = - [T d <R)] - 1 [t c <R> ; T d (R)] /* { F f w} ds (17*) 


Equation (172) gives the response state vector of a hingeless rotor blade due to a 
simple harmonic load. The hub constraint reactions, | ZTO) \ , which appear in equation 
(172) can be obtained from equation (174). It may be noted that the transfer matrix and 
all state variables are complex due to inclusion of the linear damping terms in the 
system equations of motion. 

The transfer matrix, by definition, is independent of the boundary conditions. 
Having determined the transfer matrix, the boundary conditions can be simply imposed 
to determine the response state vector as shown above for the case of a hingeless blade. 
So all the problems corresponding to various boundary conditions can be analyzed easily. 
Another advantage of the method is that it provides for the continuous integration of 
the equations. In the present study the Runge-Kutta procedure is utilized to solve the 
differential equations governing the elements of the transfer matrix (eq. (135)). One 
distinguishing advantage of the transfer matrix is that the hub constraint reactions are 
obtained as a direct result (eq. (174)). The extended transfer matrix approach is valid 
for any linear damping since normal modes are not considered. The number of stations 
at which the transfer matrix is determined depends only on the step size of the 
integration. The incremental step size can be selected arbitrarily and then the result 
compared with that obtained for a smaller increment. In this way an optimal value for 
the step size is obtained. 


RESULTS AND DISCUSSION 


The object of this study is to establish a comprehensive analytical design procedure 
for the installation of simple pendulums on the blades of a helicopter rotor to suppress 
the force and moment reactions at the hub during forward flight conditions. This 
procedure will correspond to a process of optimization with the ultimate object of 
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minimizing the root shear and moment. The minimization of these reactions wiii be 
accomplished by a systematic variation of the pendulum parameters (both geometrical 
and inertial). In this section numerical results are presented together with a discussion 
of their significance. The numerical results presented pertain to both uniform and 
nonuniform hingeless rotor blades. The properties of these blades are given in Appendix 
D. 


Simple Flapping Pendulum 

The first part of this section compares the root reactions from this study with the 
results of Murthy (ref. 15) for an undamped analysis. This comparison is made solely to 
verify the correctness of the computer program. The remainder of this section 
contains: 1) a detailed discussion of the numerical results for a 4/Rev concentrated 
load at the tip? and, 2) the results for excitation by a harmonic variation of spanwise 
airload in forward flight. 

Analysis verification . - By dropping the damping terms in the blade and pendulum 
equations of motion the case of undamped dynamic response can be recovered. The 
undamped root reactions of a uniform rotor configuration (Appendix D) are presented in 
Table 2 together with the results of reference 15. The reactions with a properly tuned 
simple pendulum are also shown in Table 2. The excitation force for all these results is 
a concentrated out-of-plane simple harmonic load at the tip of the blade with a 
frequency equal to 4/Rev. It is apparent from the data presented in Table 2 that the 
present analysis provides results which agree with the work of reference 15 for a 
uniform rotor blade without damping. The computer code has also been verified by 
duplicating published data of normal vibration modes for hingeiess and articulated rotor 
blades (ref. 16). 

Concentrated load excitation . - The force vector for a concentrated out-of-piane 
load can be written as (see eq. (168)) 

{fJ T = Looooo -F C J (175) 

where F is the amplitude of the applied load. The concentrated harmonic force at the 
tip has the following characteristics: 

Magnitude .2224 N 

Frequency 4/REV 

Both the uniform and nonuniform blades are considered for this particular excitation. 

The root reactions of the uniform blade with a flapping simple pendulum are 
computed. It may be noted that all state variables are complex due to inclusion of the 
linear damping terms. The results for the vertical shear amplitude at the hub are 
plotted in figure 9. The pendulum weight is 66.726 N (10% of the blade total weight). 
Damping at the pendulum hinge is neglected. The effect of this damping on pendulum 
tuning is discussed later in this section. 

The tuning of the pendulum is changed by varying the spanwise location of the 
pendulum for a fixed value of the uncoupled pendulum frequency. It may be noted from 
equation (53) that the uncoupled natural frequency of the pendulum is mainly dependent 
on (1) rotational speed of the blade, (2) spanwise location of the pendulum, and (3) its 


58 



TABLE 2. - CORRELATION OF UNDAMPED 
ROOT REACTIONS 




length. The simple pendulum considered here swings in the flapping plane and hence can 
generate significant out-of-plane shear force and moment at its point of connection 
with the blade. These facts can be observed from the results of figure 9. The root 
normalized vertical shear ( the shear reaction with the pendulum attached to the blade 
divided by that value without the pendulum) is altered significantly by the pendulum. 
For instance, when the pendulum is located at 15.2% of the span, the root normalized 
vertical shear is attenuated to 2.9% for a pendulum frequency equal to 146 rad/sec 
(3.87/Rev), while it is amplified to 2579% for a pendulum frequency equal to 160 
rad/sec (4.24/Rev). The latter case shows that an improperly tuned pendulum might 
amplify the root reactions rather than attenuate them, and the system would be better 
off without the pendulum. 

From the results of figure 9, it can be observed that the pendulum of frequency 146 
rad/sec (3.87/Rev) located at 18% of the span is a proper design because it attenuates 
the hub vertical shear to 1.3%. The out-of -plane moment at the hub Is reduced to 
14.2%. This choice of pendulum also attenuates the in-plane force and moment at the 
hub to 4.2% and 22.8%, respectively. 

The effect of pendulum mass on the vertical reaction is illustrated in figure 10 as a 
function of spanwise location. The results are for a simple flapping pendulum of 
frequency 146 rad/sec. It can be observed from these results that if the weight of the 
pendulum is 66 .726 N (at x « a 18% span) then the normalized vertical shear at the hub 
is 1.3%. If the weight of the pendulum is reduced by two thirds (located at 18% span) 
then the normalized vertical shear at the hub is 39.3%. This variation of the pendulum 
weight will be referred to as the weight ratio, where in this instance the weight ratio is 
33%. This result indicates that even though the uncoupled natural frequency of the 
pendulum is independent of the mass of the pendulum, a proper choice of the mass is 
also required to generate an optimum force to attenuate the root reactions. 

It can be observed from figure 9 that the minimum normalized vertical shear for 
any pendulum frequency occurs when the pendulum is located near 18% of the span. 
These values are replotted for different uncoupled pendulum frequencies in figure 11. 
Comparison of the hub vertical shear with pendulum absorber to the hub shear with no 
absorber shows the pendulum absorber to be effective in reducing hub out-of-plane 
reaction for tunings between 3.82 and 4.00/Rev. It can be recalled that the optimum 
tuning is 3.87/Rev. The steep gradient of the hub shear with pendulum tunings beyond 
optimum indicates decreasing effectiveness of the pendulum absorber and for tunings 
above 4/Rev the pendulum absorber becomes detrimental. Figure 11 demonstrates that 
in the case of an off-tuned pendulum, it is apparent that the result can be either a 
slight-to-moderate degradation in pendulum absorber effectiveness or a drastic increase 
in hub reaction depending upon the direction of tuning variation from the optimum 
pendulum natural frequency. 

Similar results have been obtained in reference 4. In this reference a semi- 
empirical analysis, using in-flight measured flap bending moment data, was developed to 
study the pendulum effectiveness in reducing 4/Rev blade root vertical shear. The 
Boeing Vertol 347 aircraft was used for the flight test. It is a tandem-rotor helicopter 
with fully articulated 4 bladed rotors. The vertical 4/Rev pendulum absorbers were 
installed on all blades of both rotors at 16% blade radius. Comparison of the predicted 
hub motion with the pendulum absorber showed that the pendulum absorber is effective 
for tunings between 3.4 and 4.0/Rev. Test results indicated that a 3.95/Rev is an 
optimum tuning. For this optimum tuning the vertical hub shears were reduced by 80 to 
90% on both forward and aft rotors. However, the vibratory blade root chordwise 
bending was increased by 25% with the pendulum installed. 

The elastic d' Uections and structural load reactions for the present uniform blade 
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corresponding to the optimum fiapping pendulum (3.87/Rev at 18%) are shown in figure 
12 together with values obtained without the pendulum. The chordwise and fiapwise 
deflections are normalized with respect to the chord (c=0.7m). It may be recalled that 
all the deflections and forces are complex. The phase angles for the deflections are 
found to be either 0° or 180° as presented in figure 12. However, the phase angles for 
the moments and forces are not exactly 0 or 180° at a few blade stations, and 
therefore it is appropriate to plot only the magnitude of these forces and moments. 

It can be observed from figure 12 that the pendulum causes a spanwise 
redistribution of the structural loads such that the hub reactions are attenuated. Figure 
12 shows that the optimum pendulum decreases the loads inboard of the pendulum and 
increases the loads outboard of the pendulum. It can also be observed that the maxima 
of the deflections, twisting moment, and vertical shear are higher than the 
corresponding quantities without the pendulum. This is the penalty one must pay for 
achieving the drastic reductions in the root reactions. The vertical shear is 
characterized by a large change at the pendulum attachment point due to the vertical 
inertial reaction of the flapping pendulum mass. This vertical reaction, which is 
imposed at the pendulum attachment point, is the basis of the blade-mounted flapping 
pendulum concept which has the purpose of attenuating the root force and 
moment reactions. 

In general, a properly tuned pendulum can attenuate the vibratory loads in two 
ways: (1) by eliminating the resonant responses of the blade by displacing the initial 
natural frequencies that are in the vicinity of the excitation frequency, and (2) by 
generating appropriate forces at its attachment point with the blade. These forces 
redistribute the loads on the blade so that the root forces are attenuated as discussed 
above. 

Figure 13 shows the hub vertical shear for pendulum tuning of the nonuniform 
blade. The blade rotor speed and precone angle are equal to 289 rpm and 0.1 rad, 
respectively. The pendulum weight is 133.452 N (13.8% of the blade total weight). The 
results indicate that a pendulum with uncoupled natural frequency of 3.85/Rev located 
at 22% of the span is a proper design because it attenuates the hub vertical shear to 
2.2% and the out-of-plane moment to 3.4%. This choice of pendulum attenuates the 
root in-plane force to 47.6%. However the root in-plane moment is amplified by 75%. 
This adverse effect can be reduced by offsetting the pendulum hinge ( ri ^) from the 
elastic axis. The effect of the pendulum hinge offset is discussed later in this section. 

For the optimum pendulum (w = 3.85/Rev and x . = 22% span) the root forces are 
computed for different values of tne pendulum mass. The results are shown in figure 
14. It can be observed that a 50% weight ratio (66.72 N) is an optimum choice. For this 
weight the normalized vertical shear at the hub is practically zero as compared to 2.2% 
with a pendulum of weight 133.452 N. This demonstrates, once again, that a proper 
choice of pendulum mass is required to minimize the root forces. 

In the above analysis the damping at the pendulum hinge is neglected. However, it 
is necessary to study the effect of the pendulum damping on the root reactions once the 
optimum pendulum tuning is established. The hub shears and moments are computed for 
different damping ratios, and these results are presented in figure 15. It can be 
observed from these results that damping of the pendulum due to friction in the pivot 
bearing has a very small effect on the hub reactions. 

Once the optimum tuning configuration and mass effectiveness are established, a 
parametric study is conducted. The parameters to be varied include prepitch, pretwist, 
precone and pendulum hinge offsets. The parametric survey is conducted for the 
nonuniform blade with the optimum tuned pendulum ( a) = 3.85/Rev and x^ = 22% 
span). The pendulum mass is 133.452 N for this study. p 
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Figure 16 shows the effect of root collective Ditch (prepitch) on the hub reactions. 
The pretwist angle for these results is equal to -9°. The prepitch has a small effect on 
the hub vertical shear. It significantly affects the hub in-plane force and moment, 
which are both increased by increasing the root collective pitch. The out-of-plane 
moment is affected moderately by changing the collective pitch. The effect of 
collective pitch on the twisting moment (not shown in figure 16) is negligible. 

The effect of the built-in twist angle (pretwist) on the hub reactions of the biade- 
mounted flapping pendulum is shown in figure 17. The root collective pitch for these 
results is equal to 13°. Figure 17 shows that the pretwist has a small effect on the hub 
vertical shear. Unlike the effect of the collective pitch, the pretwist has a significant 
effect on the root out-of-plane moment. However, for typical rotor blades (negative 
pretwist) the root out-of -plane moment is still attenuated to less than 10%. It can be 
observed from figure 17 that the in-plane force and moment are both increased by 
decreasing the blade pretwist. The pretwist has negligible effect on the twisting 
moment. It can be concluded that the flapping pendulum effectiveness in attenuating 
the root out-of-plane moment, and the in-plane force and moment is reduced by 
decreasing the blade pretwist. 

The root reactions of the blade are computed for various precone angles, and these 
results are presented in figure 18. It can be observed that the blade precone has no 
significant influence on the hub reactions. Thus, the flapping pendulum effectiveness in 
attenuating the blade root reactions is not affected by changing the precone angle. 

The results from varying the pendulum hinge chordwise offset, ru, are shown in 
figure 19, It may be recalled that the optimum flapping pendulum is n tuned with its 
hinge located at the shear center of the blade cross section ( ru - 0, £ . a 0), It can be 
argued that this is an impractical location since the shear center is usbally inside the 
airfoil envelope. The pendulum chordwise hinge offset has a small effect on both the 
root out-of-plane force and moment as shown in figure 19. Both the in-plane force and 
moment are decreased by increasing the chordwise offset. The effect of cho»dwise 
offset on twisting moment at the hub is negligible. 

Figure 20 shows the root for ces for different values of the pendulum hinge vertical 
offset, c A . It can be concluded from figure 20 that the vertical offset, 5 a, has only a 
slight effect on the hub reactions. Although not illustrated in figure 20 its effect on 
twisting moment is negligible. 

Aerodynamic load excitation . - The dynamic excitation of the rotor blade is 
provided by a quasi-steady representation of the spanwise airload distribution associated 
with the cyclic pitch environment of forward flight. As previously discussed the 
aerodynamic loads induced by the elastic deformations are also included. The 
aerodynamic load distributions consist of lift, drag, and pitching moment. Values of the 
aerodynamic parameters used to calculate the airloads are presented in Table 3. 

As previously discussed the unsteady aerodynamic forces which occur in forward 
flight are due to time dependent variations in velocity, angle of attack, and elastic 
deflections encountered by the rotating blade. The quasi-steady representation provides 
1/Rev, 2/Rev, and 3/Rev excitation airloads. The amplitude of the airloads decreases 
as the harmonic order increases. For each of the three harmonics, the amplitude of the 
vertical force is much larger than the amplitude of the in-piane force. 

For airloads which are identical on each rotor blade, the only loads and moments 
which the blades may transmit to the rotor hub are those which are integral multiples 
(N) of the number of blades (ref. 17). The N/Rev rotor hub loads are the greatest 
contributors to helicopter vibration. The source of N/Rev hub loads are N-i, N, N+l 
rotor blade airload harmonics. Therefore one solution for reducing N/Rev vibration is 
to cancel or eliminate these harmonics of blade loading. As discussed above, the quasi- 
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TABLE 3. - VALUES OF AERODYNAMIC PARAMETERS FOR THE 
NONUNIFORM BLADE IN FORWARD FLIGHT 


Advance ratio, p f 0.3 

Shaft angle, a s > deg -5 

Cosine cyclic pitch coefficient, e ic’ deg -1.5 

Sine cyclic pitch coefficient, 0^ s , deg ♦ 3.0 

Collective pitch, 6 , deg 13 

V 

Pretwist, 8 p t , deg * * "9 

Airfoil lift-curve slope, a, per rad .... 2 tt 

Profile drag coefficient, c^ Q 0.01 

Precone angle, 3 p C » rad • • , 0,1 

Blade rotational speed, , rpm ........... 289 
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steady representation of the airloads contains only three harmonics. In the present 
analysis, the effectiveness of the flapping pendulum in reducing the blade root forces 
due to 2/Rev and 3/Rev airloads is discussed. Throughout the following, unless stated 
otherwise, all numerical results pertain to the nonuniform blade. 

The root reactions of the nonuniform blade with a flapping pendulum are computed 
for a 2/Rev spanwise airload distribution. The results for the vertical shear at the hub 
are plotted in figure 21. The pendulum weight is 133.452 N. From the hub reactions it 
is observed that the pendulum of uncoupled natural frequency 2.02/Rev located at 22% 
of the span is a proper design because it attenuates the hub vertical shear to 22% and 
the out-of-plane moment to 31%, This choice of pendulum attenuates the in-plane 
moment by only 4% and amplifies the in-plane force by 66%. However for the 2/Rev 
airload the in-plane hub force is small compared to the hub vertical shear. 

For the optimum pendulum tuning ( ui s 2.02/Rev and x A = 22% span) the root 
reactions are computed for different values^of the pendulum mass and the results are 
presented in figure 22. It can be observed that a 94.53 N pendulum weight (0,71 weight 
ratio) is a proper choice because it attenuates the hub shear to 2.9 %. It also attenuates 
the hub out-of-plane moment to 17%. 

To establish the importance of representing the excitation as a distributed loading 
the responses due to the distributed airload are compared with corresponding responses 
associated with a concentrated harmonic force applied at the tip with the same 
frequency. Figure 23 shows the results for a 2/Rev concentrated out-of-plane force at 


the tip. It can be noticed from figure 23 that the flapping pendulum effect on the hub 
vertical shear is similar to the results obtained for the 2/Rev spanwise airload 
distribution (fig. 21). For specific values of the pendulum frequency and spanwise 
location, the pendulum attenuates or amplifies the root vertical shear the same way for 


both the distributed airload and the concentrated load. It is observed that an optimum 
tuning occurs at 22% spanwise location with pendulum frequency equal to 2.02/Rev. 
This optimum tuning attenuates the vertical shear and the out-of-plane moment to 17 


and 38%, respectively. It amplifies the in-piane force by 63%. It may be recalled that 
the optimum pendulum for the 2/Rev distributed airloads also occurs at x A = 22% span, 
0 ) = 2.02/Rev, and W = 94.53 N. A 

p The pendulum mass optimization for the 2/Rev concentrated load is shown in figure 
24. It can be observed that a 100.09 N pendulum weight (0.75 weight ratio) is a proper 
choice because It attenuates the hub shear to 8%. It also attenuates the root out-of- 


plane moment to 26%. 

For the 2/Rev concentrated load at the tip, discussed above, the elastic- 
aerodynamic coupling is not included. However, it is important to study the effect of 
these terms when the blade is excited by a concentrated load. Figure 25 shows the 
results for a 2/Rev concentrated load at the tip with the aerodynamic coupling terms 
included, as well as the results without these terms (W = 133.45 N). It can be noted 
from figure 25 that for both cases the root vertical shear has the same characteristics. 
However, the amount of attenuation is different. Figure 25 shows that, for the case 
with aerodynamic coupling the optimum tuning occurs at x A = 26% span and 
(jj s 2.02/Rev. This choice of pendulum attenuates the root out-of-plane shear and 
moment to 38% and 84%, respectively. The results of varying the pendulum mass 
indicate that the reaction attenuation varies only by 1.0%, The results also show that a 
133.45 N is a proper choice. Comparing the pendulum tuning for the 2/Rev distributed 
airloads, the 2/Rev concentrated load, and the 2/Rev concentrated load with 
aerodynamic coupling, it can be concluded that a pendulum tuned for a concentrated 
simple harmonic load acting at the blade tip can serve the purpose of a distributed 
simple harmonic loading at the same frequency as far as the pendulum tuning and mass 
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optimization is concerned. However, it is necessary to utilize appropriate distributed 
airloads to accurately determine the attenuation of the root reactions. 

Figure 26 shows the effect of pendulum tuning on the root vertical shear for 3/Rev 
spanwise airload distribution. The results of figure 26 indicate that a flapping pendulum 
of frequency 3.10/Rev located at 26% of the span is an optimum tuning. This choice of 
pendulum attenuates the root vertical shear and the out-of-plane moment to 4.2 and 
26.8%, respectively. It attenuates the root in-piane moment only to 97%. However, the 
in-plane force at the hub is amplified by 44%. This adverse effect is the penalty being 
paid for achieving the drastic reductions In the root out-of-plane force and moment. 

The root reactions are computed for various masses of the optimum pendulum. 
Results for the vertical shear are presented in figure 27, These data show that a 155,69 
N pendulum weight (1.17 weight ratio) is a proper choice. For this weight the vertical 
shear and out-of-plane moment are attenuated to 3.3 and 25.2%, respectively. 

The above results for a simple flapping pendulum on the nonuniform blade are 
summarized in Table 4. 


Simple Lead-Lag Pendulum 

As discussed in the previous section, the simple flapping pendulum can be tuned to 
achieve significant reductions in the root out-of-plane forces. To significantly 
attenuate the in-plane forces a simple lead-lag pendulum is needed. In this section, the 
effect of mounting a lead-lag pendulum on a rotor blade is discussed. The excitation 
loads to be considered are 4/Rev concentrated load at the blade tip and 2/Rev and 
3/Rev spanwise distributed airloads, Throughout the following all numerical results 
pertain to the nonuniform blade. 

Concentrated load excitation . - The root reactions of the nonuniform blade with a 
simple lead-lag pendulum are computed. The results for the in-plane force amplitude at 
the hub are plotted in figure 28. The study of the effect of the pendulum tuning on the 
root in-plane reaction is of prime importance since the lead-lag pendulum attenuates 
mainly the in-plane forces. The pendulum weight is equal to 133.452 N. From the 
results it can be observed that the pendulum of uncoupled natural frequency 3.85/Rev 
located at 26% of the span is a proper design because it attenuates the hub in-plane 
force to 5% and the in-plane moment to 44.6%. This choice of pendulum also 
attenuates the root vertical shear and the out-of-piane moment to 20.6 and 16.8%, 
respectively. In principle, it is feasible to design a lead-lag simple pendulum to achieve 
significant reductions in the hub in-plane force and moment without producing an 
adverse effect on the out-of-plane forces. 

For the optimum lead-lag pendulum (to * 3.85/Rev and x. - 26% span) the root 
forces are computed for different values of me pendulum mass. The results are shown 
in figure 29. It can be observed from figure 29 that if the pendulum weight ratio is 
reduced from 1.0 to 0.3, then the normalized in-plane force at the hub increases from 
5% to 39%. This result indicates that even though the uncoupled natural frequency of 
the lead-lag pendulum is independent of the mass of the pendulum, a proper choice of 
the mass is required to minimize the root forces. It may be recalled that this 
conclusion was also observed for the simple flapping pendulum. 

It is appropriate, at this stage, to compare the results for the optimum tuned 
flapping and lead-lag pendulums. Shown in Table 5 are the optimum pendulum 
parameters for both the flap and lead-lag pendulum configurations. Also shown in this 
table are the normalized root forces of the nonuniform blade due to a 4/Rev 
concentrated load at the tip. From Table 5 it can be observed that for both the flap and 
lead-lag pendulums the uncoupled frequencies are the same. Also the uncoupled 
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frequency (3.85/Rev) is in the vicinity of the external load frequency (4/Rev). It may be 
concluded that, in general, the optimum pendulum frequency depends mainly on the 
excitation frequency. This observation is due to the fact that the optimum tuned 
pendulum attenuates the hub reactions by eliminating the resonant responses of the 
blade (in addition to generating an appropriate reaction force at the hinge). The 
optimum pendulum eliminates these resonant responses by displacing the natural 
frequencies that are in the vicinity of the excitation frequency, and since the pendulum 
is more effective in altering the blade frequencies adjacent to the uncoupled pendulum 
frequency (ref. 7), then the uncoupled pendulum frequency must be close to the 
excitation frequency. As expected, the flapping pendulum significantly attenuates the 
hub vertical shear, whereas the lead-lag pendulum reduces the hub in-plane force 
drastically as shown in Table 5. Finally, it can be observed from Table 5 that for this 
particular blade the flapping pendulum produces an adverse effect on the twisting and 
in-plane moments, whereas the lead-lag pendulum attenuates the twisting and out-of- 
plane moments. 

The spanwise variations of deflections and forces for the nonuniform blade 
corresponding to the optimum lead-lag pendulum, which is located at 26% span, are 
shown in figure 30 together with the values obtained without the pendulum. The lead- 
lag and flapwise deflections are normalized with respect to the chord (c = 0.53 m). Only 
the magnitude of the forces and moments are plotted. It can be observed from figure 
30 that the lead-lag pendulum redistributes the structural loads of the blade such that 
the hub f t rees are attenuated. The pendulum decreases the loads inboard of its location 
and increases the lv«.ds outboard. Figure 30 shows that the lead-lag deflection of the 
blade-pendulum system is significantly reduced. Also, it can be observed that the in- 
plane force, V , is characterized by a large change at the lead-lag pendulum attachment 
point due to tne intertial reaction of the lead-lag pendulum mass. It may be recalled 
that the flapping pendulum produces a large change in the vertical shear. These 
reactions at the attachment point redistribute the loads on the blades in such a way that 
only a smail portion of the shears and moments are transmitted to the hub. 

A parametric study is conducted for the optimum lead-lag pendulum 
(w = 3.85/Rev and x « = 26% span). The parameters varied include prepitch, pretwist, 
precone, and pendulum ninge offsets. 

Figure 31 shows the effect of the root collective pi£ch (prepitch) on the hub 
reactions. The pretwist angle for these results is equal to -9°. It can be observed from 
figure 31 that increasing the prepitch decreases the hub in-plane force to a minimum 
value after which it increases. The in-plane hub moment is increased by increasing the 
prepitch (in the normal range). Figure 31 shows that the prepitch has a negligible effect 
on both the out-of-plane shear and moment at the hub. 

The effect of pretwist on the hub reactions is shown in figure 32, The root 
collective pitch for these results is equal to 13°. Figure 32 shows that both the hub in- 
plane force and moment are increased by decreasing the blade pretwist. Thus, the lead- 
lag pendulum effectiveness in attenuating the root in-plane force and moment is 
reduced by decreasing the blade pretwist. It can be observed from figure 32 that the 
effect of the pretwist on both the out-of-plane force and moment at the hub is small. 

The root reactions as a function of the blade precone angle are computed, and the 
results are presented in figure 33. It can be observed that the precone angle has no 
significant influence on the hub reactions. It may be recalled that the same conclusion 
was observed for the flapping pendulum. 

The results from varying the lead-lag pendulum hinge chordwise offset (n J are 
shown in figure 34. It may be noted that the optimum lead-lag pendulum is tunecT with 
its hinge located at the shear center of the blade cross section ( n ^ = 0, = 0). It can 
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be observed from figure 34 that the pendulum hinge chordwise offset has only a slight 
effect on both the hub in-plane force and moment. However, both the hub out-of-plane 
force and moment are decreased slgificantly to practically zero, after which they are 
increased, by increasing the chordwise offset. 

Figure 35 shows the root reactions for different values of the lead-lag pendulum 
hinge vertical offset, C*. This offset has a small effp • the root in-plane moment. 
However, the in-plane force at the hub is increased by i-wi* nsing the vertical offset as 
shown in figure 35. Also, both the root out-of-plane force aid moment are increased by 
increasing the lead-lag pendulum hinge offset, s *. It can also be noted that both the 
chordwise and vertical offsets have a negligible effect on the hub twisting moment 
reaction. 

Aerodynamic load excitation . - The root reactions of the nonuniform blade with a 
lead-lag pendulum are computed for a 3/Rev spanwise airload distribution. The results 
for the root in-plane force are shown in figure 36. It can be observed from figure 36 
that for this particular load the lead-lag pendulum attenuates the root in-plane force 
only by a small amount. However, the root in-plane force is small for all three 
harmonics of the spanwise airload distribution. The lead-lag pendulum of uncoupled 
frequency 3.10/Rev located at 19% of the span is a proper choice because it attenuates 
the hub in-plane force to 65%. However, it amplifies the other root reactions by 5 to 
25%. 

The effect of lead-lag pendulum tuning on the root in-plane force for a 2/Rev 
spanwise airload distribution and a 2/Rev concentrated out-of-plane load is shown in 
figures 37 and 38. It can be observed from these figures that for specific values of the 
pendulum frequency and spanwise location, the pendulum attenuates or amplifies the 
root in-plane force the same way for both the distributed airload and the concentrated 
load. The pendulum of frequency 2.01/Rev located at 9% span attenuates the root in- 
plane force to 60% for the distributed airload and 42% for the concentrated load at the 
tip. As concluded earlier, for the flapping pendulum, a pendulum can be tuned for a 
concentrated simple harmonic load at the tip as well as by a distributed simple 
harmonic loading of the same frequency. 

CONCLUDING REMARKS 

Comprehensive analytical design procedures for the installation of simple 
pendulums on the blades of a helicopter rotor to suppress the root reactions were 
developed. To achieve this goal, a frequency response analysis of typical rotor blades 
excited by a harmonic variation of spanwise airload distributions as well as a 
concentrated load at the tip was conducted. A single nonuniform rotor blade with 
hingeiess hub restraint undergoing coupled flapwise bending, chordwise bending, and 
torsional vibrations was considered. The equations of motion of both flap and lead-lag 
pendulums were derived. Expressions for the aerodynamic loads associated with the 
elastic deformations and the cyclic pitch environment of forward flight were presented. 
The transfer matrix method was used to determine the optimum pendulum tuning and 
mass to suppress the hub reactions. Lastly, the effects of various structural dynamic 
properties of the blade on the optimum pendulum configuration were established. 

The results of the study indicate that the following conclusions can be drawn. 

(1) In general, it is feasible to design a simple pendulum to suppress the vibratory 
loads at the hub of a rotor blade. Further: 

(i) A properly designed flapping plane simple pendulum attenuates the root out- 
of-plane force and moment significantly. In principle, the reductions in the root out-of- 
plane reactions can be achieved without producing adverse effects on the root in-plane 
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forces. However, in some of the cases investigated, the hub in-piane forces are 
amplified with the flapping pendulum. This is a design penalty which results from 
reducing the out-of-plane reactions. 

(ii) An optimum lead-lag pendulum design attenuates the root in-plane force and 
moment. 

(iii) For the cases investigated, the flapping plane simple pendulum is more 
effective in attenuating the root out-of-plane reactions than the lead-lag pendulum in 
suppressing the in-plane forces at the hub, 

(2) A properly tuned pendulum can attenuate the vibratory loads in two ways: (1) 
by eliminating the resonant responses of the blade by displacing the initial frequencies 
that are in the vicinity of the excitation frequency, and (2) by generating appropriate 
forces at its attachment point with the blade. These forces redistribute the loads on 
the blade so that only a small portion of the shears and moments are transmitted to the 
hub. 

(3) For optimum tuning, the parameters to be determined are the pendulum 
uncoupled natural frequency, the pendulum spanwise location, and its mass. Further: 

(i) The pendulum uncoupled natural frequency is a function of the (a) rotational 
speed of the blade, (b) spanwise location of the pendulum, (c) its length, (d) steady-state 
deflection of the pendulum, and (e) blade pitch angle at the pendulum hinge location. 

(ii) The optimum pendulum frequency is in the vicinity of the excitation 
frequency. In general, both the flapping plane pendulum and the lead-lag pendulum have 
the same optimum frequency for the same excitation frequency. In the case of an off- 
tuned pendulum the results range from a slight-to-moderate degradation in pendulum 
absorber effectiveness to a drastic increase in hub vibration. 

(iii) The uncoupled natural frequency of the pendulum is independent of the 
mass of the pendulum, and a proper choice of the mass is required to generate an 
optimum force to minimise the root reactions. 

(4) A pendulum can be tuned and its optimum mass determined by excitation with a 
concentrated simple harmonic load at the tip. However, it is necessary to utilize 
appropriate distributed airloads to accurately determine the attenuation of the root 
reactions. 

( 5 ) Damping at the pendulum hinge has a very small effect on the hub reactions 
once the optimum pendulum tuning is established. 

(6) The parametric study for a 4/Rev flapping plane pendulum and a *f/Rev lead-lag 
pendulum indicates: 

(i) Blade precone has no significant effect on the hub reactions for both types of 
pendulum. 

(ii) The effect of collective pitch on the hub vertical shear is small, but has an 
appreciable effect on the hub out-of-plane moment for a flapping plane pendulum. For 
the lead-lag pendulum the collective pitch has a negligible effect on both the hub out- 
of-plane force and moment. Increasing the collective pitch, the hub in-plane force and 
moment are significantly raised in the case of the flapping pendulum, whereas they are 
lowered to a minimum value after which they are increased in the case of a lead-lag 
pendulum. 

(iii) The pendulum effectiveness in attenuating the hub in-plane force and 
moment is reduced by decreasing the blade pretwist for both types of pendulum. In the 
case of the flapping pendulum, pretwist has a significant effect on the hub out-of-plane 
moment, 

(iv) The pendulum hinge chordwise offset has a small effect on the hub out-of- 
plane reactions in the case of a flapping plane pendulum. Also, for this type of 
pendulum the hub in-plane reactions are slightly decreased by increasing the chordwise 
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offset. In the case of a lead-lag pendulum, the hinge chordwise offset has only a slight 
effect on the root in-plane reactions, whereas the out-of-piane forces are significantly 
decreased by increasing the chordwise offset. 

(v) The pendulum hinge offset normal to the chord has only a slight effect on 
the hub reactions for the flapping piane pendulum. In the case of a lead-lag pendulum, 
the hub in-plane force and both the out-of-plane force and moment are increased by 
increasing the hinge offset. But the hinge offset has only a small effect on the hub in- 
piane moment. 

(vi) For both the flapping and lead-lag pendulum configurations the pendulum 
hinge offset in both the chordwise and normal directions had a negligible effect on the 
twisting moment reaction at the hub. 


APPENDIX A 

COORDINATE SYSTEMS AND TRANSFORMATIONS 

Several coordinate systems are used in deriving the blade and pendulum equations 
of motion. Those which are common to both the dynamic and aerodynamic aspects of 
the derivation are shown in figures A-l and A-2. The orthogonal axes system X, Y, Z 
(fig. A-l) and associated unit vectors T, 3 , K are fixed in an inertial frame &0 • It is 
assumed that the rotor hub has neither translational nor angular motion. Orthogonal 
axes x , y, Z are fixed in a inference frame 6 $ which rotates with respect to £%.at 
constant angular velocity flK. Point O, which is common to both the inertial and 
rotating frames, is located at the root of the beam. The x axis, which lies along the 
elastic axis of the undeformed beam, is inclined to the plane of rotation (and to the x 
axis) by ^th| ^recone angle 3 . The orthogonal axes x, y, z and the corresponding uni 1 ? 

vectors 1 , j , k, therefore, are*also fixed in . Beam bending deformations shown in 
figure A-2 for 0=0 are described by the displacements of the elastic axis u, v, w 
parallel tot, % k respectively. 

The transformation between the X, Y, Z coordinates of the inertial reference 
frame^^ and the x , y, Z coordinates of the rotating reference frame, (jo, in given by 
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Rotation of the x , y, Z coordinate system by the precone angle 3 about the y axis 
yields the x, y, z ?indeformed coordinate system. These two systems^are related by the 
following transformation. 
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For small values of the precone angle equation (A-2) can be approximated as 
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The relationship between the deformed x', y’, z' coordinates and the undeformed x, 
y, z coordinates illustrated in figure A-2 has been developed in refer|n^e 10. A 
transformation matrix, [ T ] , which relates their respective unit vectors i’, j',k' and i, 
ft can be defined as 
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where the transformation matrix to first order (second and 
discarded) is given by 
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higher order terms are 
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APPENDIX B 


ACCELERATION COMPONENTS 


The purpose of this appendix is to formulate the components of absolute 
acceleration as measured in the deformed coordinate system for a general point 
(x i , y,, z,) in the undeformed system, It is desired to use the deformed coordinate 
system since it is easier to formulate the pendulum equilibrium equations in this system. 
These acceleration components are to be described in terms of section coordinates and 
the elastic blade displacements and rotations to the first order. 

A position vector in the inertial reference system can be written as 

'r = xl + YJ + Zft (B-l) 


and the corresponding acceleration components are 

a^ = X ; ay s Y j a^ = Z 


(B-2) 


The components of acceleration in the rotating, nonpreconed system can be obtained 
from the component transformations in equation (A-i) as 
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X cos ft t + Y sin ft t 

Y cos ftt - X sin ft t 
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(B-3) 


Substitution of equation (B-3) into equation (A-3) yields the acceleration components in 
the undeformed system which is rotating at a small precone angle as illustrated in 
figure A-l. These components become 

** «• M 

a x = X cos ftt + Y sin ftt + Z 0p C 

• » •* 

a y = Y cos ftt - X sin ft t 

a = Z - (X cos ft t + Y sin ftt) 0 

4 Pw 

The position vector of equation (B-l) can be written for a general point (x., y,, 
Zj) in the undeformed system by successive applications of equations (A-l) and (A-3) as 

r = (xj cos ftt -y^ sin ftt -z^ 0 pc cos ftt)t 



+(Xj sin ftt +y 1 cos ftt -z. 0 pc sin &t)3 
+(x l 0 pc + (B " 5) 
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(B-6) 


The second time derivative of the components of equation (B-5) are 
X = Xj cos fit -yj sin fit -z A 3 pc cos fit 

-2 fi x^ sin fit -2 fiy^ cos fit + 2 fi Zj 0 pc sin fi t 

- fi 2 x. cos fit fi 2 y. sin fit + fi 2 z, 0 „ cosfi t 

1 l pc 

Y = x^ sin fit + y^ cos fi t -Zj 0 pc sin fit (B-7) 

+2 fi Xj cos fi t -2 fiy^ sin fi t -2 fi z^ 0 pQ cos fi t 

2.2 2 

- fi Xj sm fi t - fi yj cos fi t + fi z^ 3 pc sin fi t 

Z = 0 QC + \ (B-8) 

Substituting these derivatives into the acceieration components of equation (B-4) and 
negiecting second degree terms in the precone angle yields 

a x = x i “ fl 2 x i + z i 0 pc | 

a y = y 1 +2fix 1 -2 fi 0 pc - fi 2 y { \ (B-9) 

a z n \ +2 0^ 0 pc + fi 2 Xjl 0 pc 1 

In the deformed body, the coordinates x,, y,, z. of the general point in the 
undeformed coordinate system can be expressed in terms of the elastic deflections u, v, 
w, <j> and the section coordinates £, n (same as x', y', z’) with the aid of equation (A- 
4) as 
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These coordinates become 

Xj a X + U + £ -V' £n COS (6 + <{> ) Sin (0 + (J) ) 
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yj = v + v 1 £ + r f cos (0 + (J>) -5 sin (0 + <j) ) (B-12) 

z A s w + w' C + n sin (0 +<{>)+ 5 cos (0 + <J>) (B-13) 

When these coordinates and their time derivatives are substituted into equation (B-9) 
the acceleration components to first order become 

a = | -2 ft v - ft 2 ( x + K ) 

A 

+2 ft 5 sin (0 + <j) ) -2 ft n cos ( 0 + § ) (B-14 ) 

a y = v+ 2 ft£-ft 2 v 

+ ( ft 2 s - 5 ) sin ( e + <|> ) - ( ft 2 n - n ) cos (0 + ^ ) 

(B-15) 

a -w + ft 2 3 _ ( x + C ) 

+ n sin ( 0 + <t> ) + \ cos ( 0 + (j> ) ( B- 16 ) 


In the above acceleration expressions, the coordinates 5, n ,5 are of a general point off 
the blade. These acceleration components in the undeformed system can be 
transformed into the deformed system components by again applying the transformation 
of equation (A -4) as 
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The final acceleration components in the deformed system can be further simplified by 
applying the trigonometric approximations of equation ( 8 ) to yield 
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a y , s ri + (fi 2 xv+w + ft 2 0 pc x) sin 0 
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APPENDIX C 

AERODYNAMIC COUPLING COEFFICIENTS 

In equation (132) the elastic defer mations-dependent aerodynamic loads are 
expressed in a matrix form* The nonzero elements of these matrices are 
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For n = 1: 
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APPENDIX D 


ROTOR CONFIGURATIONS 

Two rotor configurations have been considered in this investigation. The first is a 
uniform blade, and the second is typical of an operational blade with nonuniform 
properties. 

The uniform blade has the following properties: 

Radius, in 260 

Linear twist, deg -10 

Weight per unit length, lb/in 0 . 5796 

Flapwise mass moment of inertia per unit length, lb-in 0 . 3476 

Chordwise mass moment of inertia per unit length, lb-in ....... 15 . 456 

2 8 

Flapwise bending stiffness, lb-in 0.3x10 

2 10 

Chordwise bending stiffness lb-in .0.1x10 

9 8 

Torsional rigidity, lb-in" . 0.2x10 

Mass axis of f set, in . -0.6 

Rotational speed, rpm 360 

The nonuniform blade has the properties listed in the following table and 
illustrated in figures D-l through D-6. 

Radius, in 288 

Root cutout, in 39 

Precone, rad 0.1 

Rotational speed, rpm 289 
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Figure D-2. - Rotor blade twist. 






Figure D-4. - Rotor blade chordwise c.g. distribution. 
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Figure D-6. - Rotor blade lead-lag moment of inertia. 
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Figure 2. - Equilibrium forces and moments. 






Figure 5 . - Free-body diagram of the element containing flapping pendulum. 
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Figure 6. - Lead-Jag pendulum after 4 deformation. 
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Figure 7. - Rotor blade airfoil inflow geometry. 
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Pendulum spanwise location / span 

Figure 9. - Root vertical shear for 4/Rev concentrated excitation, 
uniform blade, W a 66.73 N, n A = ? A = 0. 
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Figure 10 - Root vertical shear for 4/Rev concentrated excitation, 
uniform blade, W , = 66.73 N, h A = C * = 0. 
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(a) Lead-lag deflection. 



Radial blade station/span 
(b) Flapwise deflection. 


Figure 12. - Deflections and loads for 4/Rev concentrated excitation 
uniform blade, x A /R = .18, w = 3.87/Rev. 
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Damping ratio, i 


Figure 15. - Pendulum damping effect for <f/Rev concentrated 
excitation, nonuniform biade. 


Normalized root in-plane shear, V , and 
in-plane moment, M . y 





Frepitch, $ Q , deg 


Figure 16 ~ Root collective pitch effect for 4/Rev concentrated 
excitation, nonuniform blade. 
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Figure 17. - Pretwist effect for 4/Rev concentrated 
excitation, nonuniform blade. 
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Figure 20. - Pendulum vertical hinge offset, C «, > effect for 4/Rev 
concentrated excitation, nonuniform blade. 
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Normalized root vertical shear 
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Figure 21. - Root vertical shear for 2/Rev distributed airload excitation, 
nonuniform blade, W = 133.45 N, n A « = 0. 
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2 22. - Root vertical shear for 2/Rev distributed 
airload excitation, nonuniform blade, 

W„ 133.45 N, 0 A = C A »0. 
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Figure 23. - Root vertical shear for 2/Rev concentrated excitation, 
nonuniform blade, W = 133.45 N, ru s r . = 0. 
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Figure 24. - ^‘vertical shear for 2/Rev concentrated 
,* cltatl °n, nonuniform blade, 

W ref = 133.45 N, n 
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Figure 25. - Root vertical shear for 2/Rev concentrated excitation, 
with and without aerodynamic coupling, nonuniform 
blade, W » 133.45 N, n A = C A * 0. 
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Figure 26. - Root vertical shear for 3/Rev distributed airload excitation, 
nonuniform blade, W * 133.45 N, = 0. 








Figure 27. - Root vertical shear for 3/ Rev distributed airload 
excitation, nonuniform blade, 

W rei = .133.45 N,n A =C A = °- 
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Without lead-lag pendulum 
With iead-lag pendulur^u.. 


(a) Lead-lag deflection. 


Radial blade station/span 
(b) Flapwise deflection . 

Figure 30. - Deflections and loads for 4/Rev concentrated excitation, 
nonuniform blade, x A /R = .26, w s 3.85/Rev. 
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Figure 38. 


Pendulum spanwise location/span 

Root in-plane shear for 2/Rev concentrated excitation, 
nonuniform blade, W = 133.45 N, n A = C n = 





